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Experiments are described on the spreading upon water of drops of oil containing small 
amounts of alcohols, esters, ketones, phenols, non-ionized acids, and non-ionized amines. It is 
shown that such polar molecules do not cause spreading unless the molal concentration exceeds 
a certain value, ¢o, which varies with the nature of the polar group, the chain length, the 
possible geometrical configurations of the molecule, and with the solubility in oil. The cause of 
their spreading is demonstrated to be a reversible adsorption at the oil-water interface asso- 
ciated with the osmotic pressure of the polar molecules dissolved in the oil. The average lifetime 
of adsorption of such molecules is shown to be only a small fraction of that of ionized molecules. 
Ao, the interfacial area per adsorbed molecule at the film pressure of 11.3 dynes per cm, is 
given for a considerable number of saturated and unsaturated straight- and branched-chain 
acids and alcohols. The observed changes of Ao with chain length and polar group are explained 
on the basis of structural and kinetic considerations. Certain new spreading phenomena due to 
the formation of colloidal suspensions and crystalline precipitates of the polar compounds are 
described, and use is made of the latter to measure the solubility in oils of the various straight 
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chain saturated acids and alcohols. 





1. INTRODUCTION 


T was demonstrated in Part I! that open-chain 
fatty acids and amines having over 13 carbon 
atoms per molecule adsorbed completely at the 
oil-water interface if the pH of the water was 
adjusted to values greater than 10.5 or less than 
3, respectively. Stated in another way, the life- 
time of adsorption of the ionized polar molecule 
was found to be infinite. 

The investigations here described were con- 
cerned with the spreading of oil drops containing 
non-ionized polar molecules. In order to obtain a 
theoretically acceptable interpretation, the al- 
cohols were first observed, since in the pH range 


*W. A. Zisman, J. Chem. Phys. 9, 534 (1941). 


used (pH 2 to pH 13) no ionization could occur 
at the interface. This was followed by .the study 
of the spreading due to fatty acids on acid water 
and fatty amines on alkaline water with the pH 
adjusted in each case practically to stop ioniza- 
tion at the oil-water interface. Finally, the other 
non-ionizable molecules were studied. 

The white mineral oil used throughout the 
experiments of Part I was also employed here. 


2. SPREADING DUE TO MONOHYDROXY ALCOHOLS 


Each of the members of the open-chain satu- 
rated alcohols from hexanol (C¢) to eicosanol 
(Coo) was dissolved in the mineral oil and drops 
of each solution were spread on the clean surface 
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of distilled water at 25°C. It was found that 
spreading did not occur unless the weight con- 
centration, w, exceeded a certain value wo. The 
greater the difference w—wo, the more rapid the 
rate of spreading and the greater the maximum 
area attained. For large values the films became 
thin enough to show uniform interference colors, 
in which case they broke up soon after attaining 
their maximum areas. The lower the molecular 
weight of the alcohol, the more rapid the rate of 
spreading, and the more violent the break-up of 
the oil disk. When the maximum spreading area 
per gram of oil (zR?/m) was plotted against 
W—Wpo, a straight line was obtained for each 
substance studied (see Fig. 1). 

These phenomena are manifestations of short- 
lifetime adsorption at the oil-water interface. 
Every molecule of alcohol which adsorbs at the 
interface does not remain there indefinitely, but, 
due to the effects of the oil solubility of the 
hydrocarbon chain of the molecule and to the 
effects of the collisions of the interfacial molecules 
with one another and with the neighboring water 
and oil molecules, it leaves after a brief stay and 
disappears into.the volume of the oil drop. At 
any instant, however, the interfacial area will 
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contain adsorbed molecules, and their thermal 
agitation and mutual forces will result in the 
creation of a two-dimensional interfacial film 
pressure which will always act to increase the 
interfacial area. However, the spreading cannot 
commence until the film pressure, F, exerted by 
the adsorbed molecules exceeds the resultant of 
the surface tension forces acting on the thin oil 
disk. This resultant, Fo, was shown in Part I, 
Section 2 to be —11.3 dynes/cm for the mineral 
oil used. 

No disk spreading can occur when w is less 
than wo because the film pressure due to the 
adsorbed molecules is too small. Once spreading 
of the oil drop can take place, the alcohol 
molecules adsorb as fast as new interfacial area 
is made available. The newly adsorbed molecules 
come from the alcohol distributed within the 
volume of the oil drop in the neighborhood of the 
interface, and the volume concentration in the 
drop decreases accordingly. Finally, the spread- 
ing ceases when the weight concentration in the 
drop has decreased to the value wo, and equi- 
librium results. As evidence for the correctness of 
this description, when a multiple dropper (see 
Part I, Section 10) was employed to exert an 
external pressure on the oil drop, the disk con- 
tracted practically completely when the piston 
oil exerted an external pressure of only 15 
dynes/cm. This is quite unlike the behavior 
under the same conditions of adsorbed films of 
ionized molecules, for they contracted only as 
much as necessary to obtain the closest packing 
of the adsorbed molecules. 

Now w—w» is the amount of alcohol which 
had to flow from the oil drop to the interface to 
permit the disk to expand to its maximum 
spreading area S. Hence, if m is the mass of the 
oil drop, the number of alcohol molecules given 
up by the solution to form the film will be: 


(w—wy)m/M(6.06 X 103), (1) 


TABLE I. Aliphatic alcohols in mineral otl 
(temperature 25°C). , 





—— 








ALCOHOL N wo IN G/G coMOLES/1000G  AoIN co 
Eicosanol 20 Too insoluble —_— —_ 
Octadecanol 18 Too insoluble — en 
Hexadecanol 16 1.94 X1073 0.802 X10 57.2 x10" 
Tetradecanol 14 1.66 X1073 0.776 X10 65.7 X10™ 
Dodecanol 12 1.39 X10-3 0.749 X1072 56.2 X10" 
Decanol 10 1.14 X1073 0.722 X1072 55.4 X10” 
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where M is the gram molecular weight of the 
alcohol. If Ao is the area per interfacial molecule, 
and R is the maximum radius attained by the 
spreading oil disk, then S=7zR?, and the number 
of molecules in the interface must be: 


aR?/Apo. (2) 
Since expressions (1) and (2) are equal; 
1R?/A y= (w—wo)m/M(6.06 X 10??). (3) 


Differentiating (3) with respect to w gives the 
relation: 


d/dw(R?/m) = Ao/M(6.06X10°). (4) 


The left-hand side of (4) is the slope, tan 6, of 
the straight line whose equation is given by (3), 
and it can be obtained by plotting 2R?/m 
against w. Hence, the useful relation is obtained : 


Ao/M= (6.06 X 10**) tan 80. (5) 


This quantitative treatment may be incorrect 
if all of the polar material in the oil is not in true 
solution. If the total weight concentration of 
polar material added to the oil is w, it may 
actually consist of a colloidal portion, w., and a 
non-colloidal portion, w,, where w=w.+w,. In 
measuring Ao by diluting this solution and so 
measuring +R?/m as a function of w, the colloidal 
part, w., May vary in a complex manner depend- 
ing on the history of the dilution process. Hence 
the quantity (w—w») need not be the amount of 
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polar material adsorbed at the oil-water interface 
at concentration w. Therefore (5) can give 
incorrect values of A» under such circumstances. 
This was found to be the case in studying solu- 
tions of tri-p-tert amylphenyl phosphate in oil. 

Relation (5) does not require a precise knowl- 
edge of the state of molecular association of the 
polar material in solution, for so long as the value 
of M used is that of the unassociated molecule, 
the value of A» obtained from a spreading oil 
drop analysis will be the area per unassociated 
molecule as adsorbed at the interface. For ex- 
ample, if hexadecanol is studied and the value of 
M used in (5) is that of a single molecule, then 
Ay is found to be 57.2 K10-'* cm’. If it is assumed 
that association occurs in pairs, M is twice as 
great and so is Ag; hence the area per chain is 
unchanged. Therefore no information is given as 
to whether or not the adsorbed molecules are 
also associated in pairs when adsorbed at the 
oil-water interface. 

It was found that whereas Ao was not tem- 
perature sensitive, wo was definitely so. Since wo 
is the concentration at which F equals the 
resultant of all the surface tension forces acting 
on the rim of a flat disk of oil, and since the 
surface tensions of oil and water are affected 
markedly by temperature changes, w» would be 
expected to be temperature sensitive. In fact the 
value of wo for the long-chain alcohols was found 
to increase roughly 5 percent for a rise of 1°C. 
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This change may be due in part to the increase 
with temperature of the solubility of the alcohol 
in oil. 

The results obtained with the aliphatic alcohols 
are given in Table I and Fig. 1. In the third 
column Wp» is expressed in grams of alcohol per 
gram of oil, while in the fourth column ¢p is its 
value when expressed in gram molecules per 
1000 grams of oil. In the fifth column is given, 
Ao, the area per adsorbed molecule. 

One regularity shown by these data has been 
brought out in Fig. 2. Upon plotting co against N, 
the number of carbon atoms per molecule, it is 
found that ¢p is a linear function of NV. 

Edge loss, or the edge diffusion phenomenon 
described in Part I, Section 4, also occurred here 
during and after the spreading of the oil disk. 
This effect became more pronounced the shorter 
the chain length of the molecule. The break-up 
effects were quite characteristic of the chain 
length and concentration of the polar molecules. 

From the relation (5) a value of A» can be 
computed ; however, the existence of edge diffu- 
sion obviously may cause the calculated value 
to be too small, the error being greater the more 
rapid the rate of edge loss, or in other words, 
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the lower the molecular weight. That the value 
of Ao is nearly correct for all values of N of 
Table I will later be demonstrated by an inde- 
pendent method of measuring Ao. In Fig. 3 will 
be found a plot of Ao against N. 

It might be thought that the lower values of 
Aj» when WN decreases beyond 14 are due to the 
increased solubility in the water of the inter- 
facially adsorbed molecules. Adams’ well-known 
rule? that monolayers of alcohols and acids at 
the air-water interface are insoluble if N exceeds 
12 is not true for monolayers at the oil-water 
interface. This is to be expected since at the 
latter interface the polar molecule is attracted by 
both the water and the oil. To prove this point 
a very simple experiment was performed using 
octanol dissolved in oil at a concentration which 
did not cause spreading but did show a high 
rate of edge diffusion. When the boundary of 
the disk was compressed externally with an 
insoluble long-chain monolayer like eicosanol, or 
with a film of talc, it stopped the edge diffusion. 
After an hour the pressure was released and the 
edge diffusion commenced with undiminished 
intensity. If the octanol molecules adsorbed at 
the oil-water interface had been appreciably 
water soluble, all of them would have been 
dissolved in that time since an hour is much 
more than the time needed for all the octanol 
molecules to diffuse from the oil to the interface, 
adsorb, and finally dissolve in the water. On the 
other hand a monolayer of octanol adsorbed at 
the air-water interface dissolved in a few minutes. 

Experiments were made with mineral oil to 
which progressively larger amounts of hexa- 
decanol were added until w exceeded the solu- 
bility limit as evidenced by the subsequent pre- 
cipitation of crystals. It was found that as long 
as the solubility limit was not exceeded uni- 
formly colored films were obtained and the graph 


TABLE II. Solubilities of alcohols in mineral oil 
(temperature 25°C). 








In Sourss’s HEAVY 
PETROLATUM 


In SONNEBORN 


ALCOHOL PETROLATUM 





Octadecanol Less than 2 X10 g/g 
Hexadecanol 2.5 X10-3 
Tetradecanol 6.1 X1073 
Dodecanol 10.2 X10-3 


Less than 2 X10 g/8 
2.6 X10-3 
4.4X10 
8.8 X10 








2N. K. Adam, The Physics and Chemistry of Surfaces 
(Clarendon Press, Oxford, 1938), second edition. 
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of «R?/m versus w was a good straight line. 
When w exceeded the value w,,=2.5X10-* gram 
per gram, the graph was still a straight line, but 
numerous colored spots appeared in the oil film 
during the spreading. Experiments with tetra- 
decanol furnished similar results when w ex- 
ceeded the value of w,=6.1X10-*, and the 
colored spots spread more rapidly relative to the 
rest of the film. A similar experiment with 
dodecanol resulted in w,=10.2X10-%, but the 
incidence of the phenomenon was not as sharply 
defined as in the case of the other two alcohols. 

Since experiments showed that precipitation of 
crystals occurred only when the hexadecanol 
concentration exceeded 2.5X10-*, it was con- 
cluded that w, which was defined by the in- 
cidence of spots was actually the solubility limit. 
The resulting solubilities are in agreement with 
the well-known fact that the solubility of long 
chain alcohols in oil increases as the number of 
carbon atoms (JV) decreases. 

The appearance of these spots is due to the 
more rapid spreading of the oil disk in the 
vicinity of an undissolved crystal or of an un- 
dissolved droplet of the alcohol which has 
migrated to the oil-water interface under the 
influence of gravity. This change in the character 
of the spreading of oil drops in passing through 
the solubility limit can be used as a technique 
to measure the solubility of any slightly soluble 
polar substance which by itself is able to spread 
on water as a monolayer. However, it cannot be 
employed with much success if the solubility 
limit (Wm) is great, since the rate of disk spread- 
ing increases with the difference w—wp, so that 
before w has reached the value wm, the speed 
of spreading will have become so great that it 
will be impossible to differentiate with the 
naked eye between the spreading in the locality 
of the uadissolved crystal and that of the sur- 
rounding film. In Table II are listed the solu- 


TABLE III. Cyclic alcohols in mineral oil 
(temperature 25°C). 








ALCOHOL N 





Benzyl alcohol 7 
Diphenyl carbinol 

Phenyl benzyl carbinol 14 
Cyclohexanol 6 
1-Ethyl cyclohexanol 8 


—— 
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bilities of the three alcohols studied expressed in 
grams of alcohol per gram of oil. These measure- 
ments are probably correct to within +3 percent. 
Such accuracy was considered sufficient for the 
purpose of this work. A graph of the solubility, 
Wm, versus N is linear for Sonneborn light 
petrolatum, but is not for Squibb’s heavy 
petrolatum. 

The solubility data of the preceding paragraph 
permit an important experimental conclusion 
with regard to the mechanism involved in the 
spreading of non-ionized molecules. If the linear 
increase in the area of the spreading disk with 
increasing w were due to the fact that the inter- 
cept wo was actually the solubility limit and 
W—wW, the excess undissolved alcohol, which in 
that way was made available to be adsorbed and 
cause spreading, then necessarily the solubility 
limit would have to be greater for hexadecanol 
than for dodecanol, since from Table I wo in- 
creases with N. However, the solubility measure- 
ments of the preceding paragraph show that the 
reverse is actually the case. Hence wo cannot be 
the solubility limit, but is less than that, and the 
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difference becomes rapidly greater as N de- 
creases (compare Tables I and II). The only 
conclusion left is that the spreading of oil disks 
containing non-ionized molecules is a mani- 
festation of interfacial adsorption from a true 
solution. 

The lack of disk spreading found at 25°C for 
eicosanol and octadecanol has a simple explana- 
tion. From Table I it is evident that wo» increases 
as N increases, while from Table II the solubility, 
Wm, decreases as WN increases. In the case of 
octadecanol or of homologous alcohols of still 
greater molecular weight, w, is less than wo; 
that is, the low solubility does not permit dis- 
solving enough alcohol for the adsorbed molecules 
to exert sufficient surface pressure to cause disk 
spreading. 

As would be expected from the fact that the 
alcohols do not dissociate in water, the spreading 
effects of alcohols in oil were found to be un- 
influenced by changes in the pH of the water or 
by the presence of metallic ions. 

The spreading effects due to cyclic alcohols 
were studied much more briefly and the results 
are given in Table III. It will be noted that cp is 
greater the greater the solubility of the alcohol 
in the oil. All these alcohols showed strong edge 
diffusion, the effect being greater the smaller 
the number of carbon atoms in the longest chain 
or branch of the molecule. In all cases the edge 
diffusion was too pronounced to permit accurate 
measurements of Ao by observation of the 
maximum spreading diameters. However, wo 
could be approximately determined by seeking 
that concentration of alcohol in the oil which 
just caused the drop to begin to spread. Due to 
their insolubility in oil interesting alcohols like 
p-hydroxy diphenyl and tripheny] carbinol could 
not be studied. 


TABLE IV-A. Saturated acids in mineral oil—temp. 25°C 
(spread on 0.01 N HCl). 











AcID N WoIN G/G co MOLES/1000G  AoIN CM? 
Hexacosanoic 26 Too insoluble -- —_ 
Eicosanoic 20 Too insoluble — -- 
Nonadecanoic 19 Too insoluble — — 
Octadecanoic 18 Too insoluble — — 
Heptadecanoic 17 2.05 X10-3 0.752 X1072 63.8 X10-16 
Hexadecanoic 16 1.71 X10-3 0.667 X1072 72.0 X10-16 
Pentadecanoic 15 1.58 X10-3 0.653 X10~2 72.6 X10716 
Tetradecanoic 14 1.37 X10-3 0.601 X10~2 70.5 X10716 
Tridecanoic 13 1.30 X10-3 0.607 X10-2 68.9 X10716 
Dodecanoic 12 1.11 X10-3 0.556 X1072 69.8 X10-16 
Undecanoic 11 0.98 X1073 0.530 X1072 57.4 X10716 
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All of the alcohols described here were East- 
man preparations except the 1-ethyl cyclohexanol 
which was kindly supplied by W. A. Mosher? of 
the Hercules Powder Company and the eicosanol 
which was prepared by Dr. Homer Adkins of the 
University of Wisconsin. 


3. SPREADING DuE TO Non-IONIZED 
MONOCARBOXYLIC ACIDS 


Each of the saturated open-chain acids de- 
scribed in Part I was also studied when spread on 
0.01 normal HCl, and it was found in each case 
that precisely the same spreading phenomena 
occurred as have been described for the various 
alcohols. In Fig. 4 will be found graphs of rR?/m 
against w for these acids, while the results suit- 
able for comparison with those for the alcohols 
will be found in Tables IV-A and IV-B. The 
eicosanoic, nonadecanoic, and _ pentadecanoic 
acids were highly purified specimens supplied by 
Dr. Roger Adams of the University of Illinois; 
the oleic and elaidic acids were freshly prepared 
in a high state of purity by Dr. J. Kass of the 
University of Minnesota. The other acids were 
Eastman preparations. 


3 W. A. Mosher, J. Am. Chem. Soc. 62, 552 (1940). 
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The preceding study of the spreading of al- 
cohols furnished the clue as to the cause of the 
non-spreading of oil drops containing low con- 
centrations of long-chain fatty acids on acid 
water (see Part I, Section 13). On acid water 
the fatty acids are not ionized, hence they adsorb 
like any other non-ionized molecules and cause 
disk spreading only when w exceeds a value wo 
of the order of magnitude of 10-*. These con- 
clusions evidently extend the earlier conclusion 
of Langmuir* * which concerned the adsorption of 
stearic acid at the oil-water interface. Therefore, 
in contrast with the behavior of ionized mole- 
cules, all non-ionized polar molecules have very 
short lifetimes of adsorption. 

A graph of co against N (see Fig. 2) is a straight 
line. Also the mean value of co for acids is 
0.613 X10-* while for alcohols it is 0.762107. 
The lower value for the acids can be considered 
either as a manifestation of the smaller solubility 
in oil of the acids, or more likely, as evidence of 
the greater energy of adhesion of the carboxyl 
group to water, than of the hydroxyl group. 

Upon comparing the value of Ao for each 





*I. Langmuir, J. Frank. Inst. 218, 143 (1934). 
°I. Langmuir, Science 84, 382 (1936). 
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alcohol and the corresponding acid, it is evident 
that it is always greater for the acid. This is not 
due to the effect of a slight dissociation of the 
fatty acid at the water interface for the same 
results were obtained using acid water varying 
in strength from 0.01 to 0.05 normal. 

The lower values of co for the cis- and trans- 
isomers oleic and elaidic acids, as well as for 
undecenoic acid, than for the corresponding 
saturated acids (see Figs. 5 and 6) are note- 
worthy. An evident interpretation is that the 
presence of double bonds increases the average 
adsorption lifetime. 

As an example of the great sensitivity of Ao 
and ¢» to the state of purity of the acid prepara- 
tion used, Eastman oleic acid gave a straight 
line graph of 7R?/m against w only for large 
values of w, and from this slope the values for 
Ao and ¢o were calculated to be 4610~'* and 
0.25X10-*, respectively. In contrast the pure 
acid prepared by Dr. Kass gave an excellent 
straight line for all values of w—wo, and Ao and Co 
were calculated to be 79.8X10-'® and 0.564 
X10-*, respectively (compare Figs. 5 and 6). 
The values of A» and Co for oleic, elaidic, and 
undecenoic acids when plotted in Figs. 2 and 3 
led to the suggestive curves shown. From the 
appearance of the graph for Eastman erucic acid 
(see Fig. 6) it was concluded that polar im- 
purities were present. 

Just as in the case of alcohols, the saturated 
fatty acids showed an increasing value of co as 
the number of carbon atoms per chain (J) in- 
creased, and similarly it was found that the 
solubility in the mineral oil decreased as N 
increased. Thus the solubility of hexacosanoic 
acid (N=26) was less than 10~‘, that of eico- 
sanoic acid (V=20) was between 0.4 10-* and 
1.010-%, that of octadecanoic acid (N=18) 
was approximately 1.010-%, and that of tetra- 
decanoic acid exceeded 3.0 10-*. 

It was found that oil drops containing octa- 


TABLE IV-B. Unsaturated acids in mineral oil—temp. 25°C 
(spread on 0.01 N HCl). 











AcID N Wo IN G/G coMOLES/1000G AoIN CM? 
Erucic 22 1.22 X1073 0.361 X1072 36.8 X107'6 
Oleic 18 1.59 X10-3 0.564 X10 79.8 X10~% 
Elaidic 18 1.53 X10-3 0.541 X10~2 72.7 X10716 
Undecenoic 11 0.650 X 10-3 0.353 X10 61.1 X10~'6 
Linolenic 18 0.560 X1073 0.202 X1072 56.5 X107% 
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decanoic acid did not spread on acid water at 
25°C. This behavior was similar to that of 
octadecanol already described (Section 2) and 
had the same cause; i.e., wo for octadecanoic 
acid would of necessity, according to Table IV-A, 
be greater than for heptadecanoic acid, and 
therefore greater than 2.010-%. Since the solu- 
bility of octadecanoic acid was less than 1.0 
X10-*, disk spreading was impossible when 
using this mineral oil. 

The results obtained with solutions in mineral 
oil of members of a homologous series of branched- 
chain acetic acid derivatives are presented in 
Table V. These were kindly supplied from 
preparations made some years ago under the 
direction of Dr. Roger Adams at the University 
of Illinois. The rapid decrease in Ao for the last 
four acids was due in greater part to the rapid 
increase in the rate of edge diffusion as the 
length of the longest branched chain (N;) de- 
creased from 15 to 8 carbon atoms. It will be 
noted that co is approximately a linearly de- 
creasing function of N, (see Fig. 2). Due to the 
presence of small amounts of polar impurities, 
possibly formed in aging, more accurate measure- 
ments of wo and A» were not possible. As in 
the case of the oleic acid, already mentioned, 
the impurities caused sufficient departure of the 
aR?/m vs. w curve from linearity for low values 
of w—wW») to make it difficult to measure accu- 
rately the intercept and slope and so obtain wy 
and Ao. 

Various acid molecules containing closed hydro- 
carbon chains were also studied and the limits of 
the oil solubility are given (Table VI). The 
chaulmoogric and hydnocarpic acids were pure, 
freshly prepared specimens obtained from Dr. 


TABLE V. Aliphatic branched chain acids in mineral oil— 
temp. 25°C (spread on 0.01 N HCl). 








; Ao X1016 
AcIp N M wo X108 co X10? cm? 





*N-Pentadecyl 
acetic 0+15+2 15 2.05 0.75 63.8 
Methyl N-tetra- 


decyl acetic 1+14+4+2 14 (3.0-3.2) 1.15 (50-62) 
Ethyl N-tridecyl 

acetic 2+13+2 13 2.9 1.1 55 
N-propyl N-dodecyl 

acetic 3+12+2 12 (3.5-3.6) 1.3 (21-32) 
N-butyl N-undecyl 

acetic 4+11+2 11 (3.9-4.2) 1.5 (15-26) 
N-heptyl N-octyl 

acetic 7+8+2 8 4.9 1.8 13 








* This is heptadecanoic acid (see Table IV-A) written for purpose of 
comparison as an acetic acid derivative. 
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Roger Adams, while the other acids were East- 
man chemicals. The low values of wo obtained 
for the chaulmoogric and hydnocarpic are to be 
attributed, at least in part, to the low oil solu- 
bilities of the hydrocarbon tails of the molecules 
which consist of a cyclopentyl group at the end 
of a saturated straight chain. Since w,, was less 
than wo it was not possible to measure Ap. 
As for the last three acids of this table, like the 
closed-chain alcohols of Table III, the edge 
diffusion effects were so pronounced as to prevent 
any accurate measurement of Ao. 

A series of arylstearic acids already described 
in the literature® were kindly made available by 
Dr. Stirton of the Division of Industrial Farm 
Products Research of the Department of Agri- 
culture. In Table VII will be found the results 
of measurements of wo and co on 0.01 normal 
hydrochloric acid. It will be noted that co does 
not vary considerably, the mean value being 
0.23 10-*. The lower values of co as compared 
with those of the single chain molecules of 
Table IV-A cannot be due to the relative oil 
solubilities of the hydrocarbon portions of the 
acid molecules, for the arylstearic acids have the 
greater oil solubilities. It is suggested that either 
the lifetime of adsorption is greater for these 
arylstearic acids, or that they exert a greater 
surface pressure when adsorbed than does a 
monolayer containing the same number of stearic 
acid molecules per unit area. 


4. SPREADING DuE TO Non-IONIZED AMINES 


Although oil drops containing concentrations 
of the order of magnitude of w= 10~ of aliphatic 
amines were found to spread readily on acid 
water and not at all on alkaline water (see Part I, 
Fig. 8), spreading also resulted on water of high 
pH if w exceeded a critical concentration «>. 
Values of wo are given in Table VIII for various 
carefully purified primary aliphatic amines fur- 
nished generously by the chemical research labo- 
ratory of Armour and Company. 

The concentration co is an increasing linear 
function of N as shown in Fig. 2. Evidently the 
behavior of the aliphatic primary amines on 
alkaline water is entirely similar to that of the 


6 A. Stirton and R. Peterson, Ind. Eng. Chem. 32, 1137 
(1940). 

















THE SPREADING OF OILS ON WATER. II 


TABLE VI. Various closed-chain acids in mineral oil—temp. 25°C (spread on 0.01 N HCl). 














AcID N wo co MOLES/1000 G Wm IN G/G 
Chaulmoogric 18 (1.0-1.5) x 1074 (0.36—0.54) x 107% 1.6XK10*<wy» 
<5.3XK10~ 
Hydnocarpic 16 5.51074 0.22 x 10-3 Wm >9.5 X 10-4 
Dibenzyl acetic 7+7+2 1.41074 0.58 x 10-3 Wm >2.7 X 10-3 
Di-n-hepty] acetic 7+7+2 3.6 10-3 1.4107 Wm >6X 10-3 
Cyclohexane carboxylic 7 (1.5-1.8) x 10-3 Wm >5X 10-3 


(1.2-1.4) x 10-* 








aliphatic alcohols on water and of the aliphatic 
acids on acid water. In fact, so long as the polar 
end of the molecule cannot ionize while at the 
oil-water interface, all aliphatic polar molecules 
appear to behave analogously as regards the 
general character of the adsorption phenomena. 

When a primary amine caused an oil drop to 
spread on alkaline water, it spread to a definite 
maximum diameter, stopped a few seconds, and 
then resumed spreading over a much larger 
area to make a thin colored film in which holes 
were formed during the final stages of spreading. 
Whereas the first stage of spreading was a 
characteristic effect due to the short lifetime of 
adsorption of undissociated amine molecules, the 
second stage was caused by a chemical reaction 
between the interfacial amine film and carbon 
dioxide dissolved in the water. This effect was 
first identified by noticing that blowing one’s 
breath on the oil disk immediately accelerated 
the second stage of spreading. It was found not 
to be due to the air pressure, and it was con- 
cluded that the effect might be due to the carbon 
dioxide in the breath. Finally, when carbon 
dioxide vapor from a piece of ‘“‘dry ice’’ was 
formed over the disk, the effect showed immedi- 
ately in a most pronounced manner. When 
Na2CO; was dissolved in the water, the secondary 
spreading effects were not accelerated, but when 
NaHCO; was dissolved in the water that type 
of spreading was always greatly accelerated. 
Therefore, it was concluded that the effect was 
caused by a reaction between the amine mole- 
cules adsorbed at the interface and the high 
HCO; ion concentration in the water surface 
resulting in the formation of a highly adsorbed 
carbonate of that amine. 

A number of secondary and tertiary aliphatic 
amines were investigated in a similar fashion and 
the results are presented in Table IX. The 
higher values of wo for these amines as compared 


to the primary amines are possibly to be associ- 
ated with their much greater oil solubilities. 
Similar experiments with various aromatic and 
naphthenic amines produced entirely similar 
results (see Table X). 

Edge diffusion was found to a remarkable 
extent in secondary and tertiary amines. There 
is no doubt that in the case of the branched 
chain molecules the total number of carbon 
atoms per molecule N does not alone control the 
rate of edge diffusion, for there is a pronounced 
difference between branched chain and single 
chain molecules having the same value of JN. 
This is in agreement with conclusions reached in 
Part I, Section 4. 


5. SPREADING DUE TO ESTERS 


Results of studying the spreading of solutions 
of esters are presented in Table XI. The effects 
were quite similar to those found for the other 
non-ionized molecules. 

The tri-p-tert amylphenyl phosphate was a 
laboratory product of the Sharples Solvents 
Corporation and the tri-p-cresyl phosphate was 
an Eastman preparation. Both appeared to dis- 
solve readily in the oil, but the very gradual 
formation of tiny holes in the otherwise stable 
colored oil films formed gave good evidence for 
concluding that some of the phosphates were 
present as a colloidal dispersion in the oil. The 


TABLE VII. Arylstearic acids in mineral oil—temp. 25°C 
(spread on 0.01 N HCl). 











MOLE 
AcID WT. wo co 

Phenylundecylic 263 0.42 x10-3 0.17 X10~ 

Phenylstearic 360 0.88 X10-3 0.24 X10-2 
Ethylphenylstearic 388 1.05 X10-3 0.27 X10 
Diisopropylphenylstearic 430 1.54 X1073 0.35 X1072 
p-Chlorophenylstearic 395 0.66 X 1073 0.17 X1072 
Xylylstearic 293 1.26 X1073 0.32 X107s 
Pseudocumylstearic 406 0.85 X1073 0.21 X1072 
Cymylstearic 416 0.97 X10-3 0.23 X1072 
Tetrahydronaphthylstearic 416 0.87 X10-3 0.21 X10: 
Xenylstearic 441 0.48 X1073 0.11 X1072 





| 
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tri-amyl citrate and the tri-amyl borate were 
semi-commercial products of the Sharples Sol- 
vents Corporation, and each showed the high 
rate of edge diffusion to be expected from mole- 
cules containing only short hydrocarbon chains. 
The other esters were Eastman preparations. 

A group of arylstearic acid esters obtained from 
Dr. Stirton® was also studied. These were methyIl- 
and phenyl-undecylate, ethylphenyl stearate, 
ethylxylyl stearate, methylphenoxyphenyl stea- 
rate, butylphenoxypheny! stearate, and butyl- 
xylyl stearate. All gave values of wo in excess of 
5X10-*, which is from 5 to 10 times as great as 
the values for the corresponding arylstearic acids 
(see Table VII). This is presumably due to the 
shorter lifetime of adsorption of the ester group 
and to the much higher oil solubilities of esters 
than acids. 


6. SPREADING DUE TO KETONES 


Oil solutions of various aliphatic and aromatic 
ketones gave results similar in their general 
nature to those already found for other non- 
ionized molecules, and the results are given in 
Table XII. The ketones, presumably due to 
their shorter lifetimes of adsorption and their 
higher solubilities in oil, had much larger values 
of co than did any of the long-chain acids or 
alcohols. 

The solubilities in oil were measured only 
for the least soluble ketones, i.e., di-undecyl 
ketone (W» = 15 10-*) and di-heptadecyl ketone 
(Wm =1.5%10-%). In the case of the latter, it 
was not possible to dissolve enough of the ketone 
in the oil to reach the value wo, and hence no 
oil spreading could be obtained. 

All the ketones studied are Eastman prepara- 
tions excepting the di-heptadecyl ketone which 
was obtained from the Chemical Laboratory of 
Armour and Company. 


TABLE VIII. Aliphatic amines in mineral oil— 
temp. 25°C (spread on 0.01 N NaOH). 











AMINE N wo IN G/G co IN MOLES/1000 G 
Octadecyl amine* 18 Too insoluble == 
Heptadecyl amine 17 2.5 X10-3 0.98 X1072 
Hexadecyl amine 16 2.2 10-3 0.91 X1072 
Tetradecyl amine 14 1.5 X10-3 0.71 X1072 
Dodecyl amine 12 1.1 X10-3 0.59 X10 
Octadeceny! amine 18 0.91 X10-3 0.34 X1072 








* Solubility, wm, is between 1.7 X1073 and 2.3 X10-3 g/g. 
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7. SPREADING EFFECTS DUE TO PHENOLS 
AND NAPHTHOLS 


Adding either phenols or naphthols to oil did 
not cause spreading unless w exceeded a critical 
value wo which is listed in Table XIII. There 
was no observed effect of varying the pH of the 
water except at high pH (say pH 11) where it 
was found an increased rate of edge diffusion 
was occurring. The presence of metallic ions 
made no change in the phenomena of spreading 
or edge loss. These observations are not un- 
expected, since phenols are well known to behave 
as very weak acids. Oil soluble phenols of 
longer chain length might prove interesting but 
have not been obtainable. 

The three chlorinated phenols listed in Table 
XIII are semi-commercial products of the 
Sharples Solvents Corporation, and as would be 
expected from their high oil solubilities the 
values of co were high. The other phenols and 
the two naphthols were Eastman preparations. 


8. THEORETICAL CONSIDERATIONS 


A qualitative explanation of these adsorption 
phenomena at the oil-water interface can be 
given in terms of the kinetic concepts originated 
by Langmuir in developing the theory of mono- 
molecular adsorption of gases on solids.’ 

Let v be the rate polar molecules strike unit 
area of the interface, oo be the actual cross- 
sectional area of an adsorbed molecule, a be the 
coefficient of accommodation, 7 be the average 
lifetime of the adsorbed polar molecules, and 
the average number of polar molecules adsorbed 
on unit area under equilibrium conditions, then: 


nm av 
—=—(1—n9). (6) 


T Go 


This relation is based on the assumption that 
no forces exist between the adsorbed molecules. 
It has already been suggested here and in Part I 
that forces did exist between the adsorbed polar 
molecules and oriented hydrocarbon molecules 
derived from the oil resulting in the formation 
of a compact mixed film at the oil-water inter- 
face. In addition, evidence will be given here for 
concluding that cohesive forces between the 


71. Langmuir, J. Am. Chem. Soc. 40, 1361 (1918). 
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TABLE IX. Aliphatic secondary and tertiary amines in mineral oil—temp. 25°C (spread on 0.01 N NaOH). 






II 








AMINE N 
Di-n-amylamine 5+5 
Tri-n-amylamine 5+5+5 
Di-octylamine 8+8 
Tri-octvlamine 8+8+8 


Wo >2.2 KX 107 g/g* 
Wo>2.6X 1073 g/g* 

20 X 1073 < wo <30 X 10% 
Wo < 13.8 107% 







we co 





co>1.4X10 
co>1.1 x10 
8.3K 107? <e9 << 12.4 107 
6o<3.0K 10 











* Pronounced edge loss. 


hydrocarbon chains of the adjacent polar mole- 
cules become the most important forces when N 
exceeds 13 carbon atoms per molecule. 

If the solution of polar molecules is dilute, 
there will apply the well-known Van’t Hoff 
relation between osmotic pressure p, molal con- 
centration c, and absolute temperature 6. If the 
density of the oil is p and R is the gas constant, 
then: 


p =(cp/1000)R8. (7) 


From the kinetic theory of ideal gases if M 
is the molecular weight: 


v=kp/(Me)?. (8) 


Even in dilute solutions relation (8) may be no 
more than a rough approximation. However, in 
dilute solutions v will certainly be proportional 
to p, and, although its variation with @ may be 
complex, it will decrease as M increases, since 
the viscous resistance offered by the oil to the 
motion of polar molecules will vary with M and 
with the shape of the molecule. 

It has been shown here that + is infinite if the 
long-chain molecule is ionized at the oil-water 
interface, and is a very small fraction of a second 
if it is not ionized. It is reasonable to assume that 
if there is cohesion between the adsorbed polar 
molecules 7 will be increased, while it will be less 
the greater the solubility’of the polar molecule 
in the oil, other things being equal. 

A quantity of importance here is A, the 
average area of interface occupied by the ad- 
sorbed polar molecules. Since by definition 
A=1/n, Eq. (6) can be written: 


A=o,/avr+ap. (9) 


Evidently the graph of A versus v is a rectangular 
hyperbola with a horizontal asymptote at A =ap. 
In these experiments » was increased either by 
increasing the concentration c of the polar 
molecules in solution, or by decreasing the chain 





































length N. Actually, for high values of v the curve 
may not be hyperbolic, since under such circum- 
stances forces between polar molecules may 
become important. 

The decrease in co with N for the saturated 
and unsaturated straight-chain acids (see Fig. 2) 
cannot be due to edge diffusion, for that would 
cause Co to increase as N decreases. The most 
likely cause is the increased film pressure (for a 
given value of surface concentration) as N de- 
creases. This is in conformity with the rule 
formulated in the study of adsorption at the 
air-water interface, to the effect that in a 
homologous series, provided mutual cohesive 
forces between the polar molecules are not 
important, as N decreases the film pressure rises 
at constant area per molecule. 

The fact that the graph of co vs. N for acids is 
below that for alcohols can be attributed to a 
greater energy of adsorption of the non-ionized 
carboxyl group than the alcohol group (or a 
greater lifetime of adsorption). Therefore, for a 
given value of NV a smaller molal concentration 
is needed to obtain the pressure Fy. The graph 
for the unsaturated acids is below that for the 
saturated acids, since, for a given value of N at 
the air-water interface, an unsaturated acid 
exerts a higher film pressure than a saturated 
acid at the same surface concentration. In other 
words, the unsaturated acids are less bound by 
cohesive or van der Waals forces between the 
hydrocarbon chains. 

The acetic acid series of Table V is particularly 


TABLE X. Other amines in mineral oil—temp. 25°C 
(spread on 0.01 N NaOH). 








AMINE N wo co 





Dicyclohexylamine 12 wo >3.6 X1073* co>2.2 1072 


Diphenylamine 6+6 wo >2.9 X1073 co>1.7 X10-2 
Triphenylamine 6+6+6 wo >3.9 X1073 co>1.6 1072 
Dibenzylamine 74+7 wo >1.9 X1073* co>0.96 X10~2 
Tribenzylamine 74+7+7 wo >6.9 X10-3* co>2.4 X10? 








* Pronounced edge loss. 
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TABLE XI. Esters in mineral oil—temp. 25°C 
(spread on 0.01 N HCl). 











EsTER N wo IN G/G co MOLES/1000 G 
Ethyl myristate 16 wo>15.0 X10 co>5.9 X10? 
Methyl stearate 19 we> 6.01073 co >2.0 X10 
Ethyl stearate 20 wo> 6.1107 co>2.0 X1072 
Tri-p-cresyl 

phosphate 74+7+4+7 1.48 X10-3 0.4 X10-2 


Tri-p-tert amyl- 
phenyl phosphate 11+11+11 1.5 X1073 <wo 0.28 X10? <co 
1.9 X1073 <0.35 X1072 


Triamy] citrate 7474641 '2X10-3 <wo 5.0 X10? <co 
<2.5 X10-3 <6.2 X10 
Triamyl] borate §54+5+5 wo >7 X10-3* co>2.6 X1072 








* Showed strong edge loss. 


interesting. From Fig. 2 it is seen that the graph 
of co vs. N is a straight line with a negative slope. 
In plotting this graph the value of N used was 
the number of carbon atoms (Jj) in the longest 
branched chain of the molecule (see column 3 of 
Table V). Upon constructing models of these 
molecules and considering various possible con- 
figurations of each, it was found that for any 
given member of this series some configurations 
were, for geometrical reasons, less likely to 
permit adsorption than others. Some configura- 
tions could be described roughly as various 
forms of the letter Y with the carboxyl group at 
the base, and others were like the inverted 
letter V with the carboxyl group hanging 
from the vertex between the open legs of the 
letter. The number of such poorly adsorbed 
configurations appeared greatest for the molecule 
> CHCOOH when JN, and Ne» were nearly 
equal, and were least when N/N2 was greatest. 
Therefore the greatest molal concentration of 
acid molecules was necessary to cause oil disk 
spreading for the molecules having N; and N,2 
most nearly equal. This was in agreement with 
the experimental results just described. 

From Fig. 3 it is seen that the graph of Ao 
vs. N has a maximum at N=14 in the case of 
the alcohols, and at N= 15 for the acids. Relation 
(9) makes it clear why the early increase of Ao 
with JN is due to the decrease in v caused by the 
lowered speed of diffusion of the polar molecules 
in oil as NV increases. The decrease in Ao after NV 
exceeds 14 is due to the increased importance of 
the cohesional forces between adjacent hydro- 
carbon chains of the polar molecules, and this 
significantly happens when N approaches values 
for which the solubility in oil becomes low. This 
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also makes the parallel behavior of the curves 
of Fig. 3 for alcohols and acids readily under- 
standable, since the mechanism essentially in- 
volves reactions between the hydrocarbon chains 
rather than the polar groups. 

The curves P,;—P:2 and Pi—P; of Fig. 3 were 
obtained by plotting A» for the unsaturated 
acids (undecylenic, oleic and elaidic). The graphi- 
cal point for the first acid coincided within the 
experimental error with the point for saturated 
acids of that chain length, while the graphical 
points for oleic and elaidic acids were at P, 
and P;. Upon noting that the ordinates corre- 
sponding to points P2, P3, and P, had descending 
magnitudes in the order of cts-, trans-, and 
normal-saturated configurations for 18-carbon 
chains, it was concluded that the cause of these 
differences was in the greater tendency of the 
saturated chains to cohere. The trans-configura- 
tion would be expected to be intermediate in 
behavior since it is more nearly a straight chain 
than the cis-configuration. These geometrical 
differences are quite evident upon inspecting 
molecular models of these substances, and their 
application to the interpretation of film phe- 
nomena at the air-water interface was indicated 
by Adam and Jessop* and later by Marsden and 
Rideal.® 

If the various types of polar compounds 
investigated here are arranged in the order of 
the values of co (as in Table XIV), they can be 
placed in three groups. The same grouping is 
obtained if values of co are arranged in their 
order of magnitude for polar molecules having 
the same value of N. It will be seen that the 
relative order of the three groups is the same as 
the relative order of the solubilities of these sub- 
stances in the oil. This demonstrates the im- 
portance of the solubility in oil as a variable to 


TABLE XII. Ketones in mineral oil (temperature 25°C). 











KETONE N Wo IN G/G co IN MOLES/1000 G 
Methylnony! ketone 10 12.2 X10-3 7.2 X10" 
Di-n-amyl ketone 10 15.1 X10-3 8.9 X10? 
Di-undecyl ketone 23 12.0 X10-3 3.6 X10~? 
Di-heptadecyl ketone 35 Too insoluble _ 
Di-pheny! ketone 13 13.4 X1073 7.3 X10~? 
Di-benzy] ketone 15 6.0 X10-3 2.9 X10 
Phenylundecy! ketone 18 41.7 X10-3 16.0 X10 








8 of K. Adam and G. Jessop, Proc. Roy. Soc. A101, 516 
(1922). 
9 J. Marsden and E. Rideal, J. Chem. Soc. 1163 (1938). 







































THE SPREADING OF 
be considered in the study of the oil spreading 
due to additive polar compounds. 

A result of considerable interest is evident 
from an inspection of Fig. 3. All the saturated 
straight chain alcohols and non-ionized acids 
have values of A» ranging from 50 to 72 square 
angstroms, whereas the same polar molecules at 
the same temperatures and film pressures would 
occupy only 20 to 25 square angstroms at the 
air-water interface. One would expect A» to be 
higher in the former circumstances due to some 
degree of film disorganization resulting from the 
thermal agitation of the hydrocarbon molecules 
of the oil. However, at film pressures as high as 
11.3 dynes/cm and in proximity to high mo- 
lecular weight hydrocarbon molecules, a much 
closer packing would be more likely. This sug- 
gests the same conclusion given in Part I, i.e., 
the adsorbed polar molecules form a fairly closely 


TABLE XIII. Phenols and naphthols in mineral oil— 
temp. 25°C (spread on 0.01 N HCl). 











PHENOL N wo IN G/G co IN MOLES/1000 G 
p-tert-amyl phenol 11 0.5 X1073 0.3 X1072 
p-cyclohexyl phenol 12 0.5 X1073 0.3 X1072 
2-chloro-4-tertamyl- 

phenol 5+6 wo< 8.810" co <4.4 X1072 
2, 6-dichloro-4-tert- 

amyl phenol 5+6 wo >12.5 X10-3 co>5.3 X1072 
2-chloro-4, 6-di-tert- 

amyl phenol 5+5+6 wo >40.0 X10-3 co>1.5 K1072 
a-Naphthol 12 1.0 X10 0.6 X10~2 
8-Naphthol 12 1.0 X10-3 0.6 X10-2 
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TABLE XIV. Effect of solubility on 
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magnitude of Co. 














MAGNITUDE 
SUBSTANCE co RANGE Group OF Co 
Chaulmoogric and 
hydnocarpic acids (0.2-0.5)10-3 I 10-3 
Arylstearic acids (0.2-0.3)10-2 


Unsaturated acids 

Saturated acids 
(Straight-chain) 

Saturated alcohols 


(0.35—-0.56) 10-2 


(0.5-0.75)10~2 
(0.72-0.80) 10-2 








II 


Straight-chain amines (0.6-1.0)10~2 
Esters (0.2—.60) 107! Ill 
Ketones (0.3-1.6)107! 107! 








packed mixed film with hydrocarbon molecules 
derived from the oil. Such a film can be con- 
sidered a two-dimensional solution. Obviously 
the proportion of polar to non-polar molecules in 
the film will vary with c, and also with the 
intensity of the adhesional forces between the 
hydrocarbon chains of the polar molecules and 
of the non-polar molecules. The correctness 
of these conclusions will be demonstrated in 
Part III. 

The cooperation of Dr. Roger Adams of the 
Department of Chemistry of the University of 
Illinois has been a source of encouragement, as 
well as the essential pure organic compounds 
already listed. A wealth of useful chemical 
information concerning hydrocarbons and polar 
additives was made available during the course 
of this work through numerous discussions with 
Dr. Parry Borgstrom. 
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Determination of Diffusion Coefficients 


B. SERIN* AND R. T. ELLICKSON 
Department of Physics, Polytechnic Institute of Brooklyn, Brooklyn, New York 


(Received July 7, 1941) 


Diffusion equations derived from Fick’s law are discussed and applied to experimental data 
in the literature on the diffusion of oxygen in copper and on reaction rates for chemical reactions 
in the solid phase. Previous limitations on the applicability of the formulae are removed by 
the use of a simple graphical method so that the theory can be applied to data to which it was 
previously inapplicable. It is pointed out that the exact solution to Fick’s law obtained is 
equivalent to the parabolic law in the early stages of the diffusion process. The activation 
energy for the diffusion of oxygen in copper is found to be 50.5 kcal. per g-atom. 





INTRODUCTION 


N an earlier paper! a set of equations was 
derived governing diffusion in solids under 
various boundary conditions. It was shown, for 
example, that if a slab of infinite extent and 
thickness L, in which the concentration of a 
diffusing material is zero at t=0, is placed in a 
region where the concentration is maintained at 
a value Cy at the boundaries of the slab, one could 
obtain the following expression for a, the fraction 
of completion of the diffusion process, 


1—a=(8/2r?) > (1/n?) exp (—n?x*Dt/L*). (1) 


n odd 


In this expression ¢ is the time and D is the 
diffusion coefficient defined by Fick’s law, 


V?C=(1/D)(dC/dt). (2) 


Similar expressions were obtained for the case 
of a sphere and a cylinder, the expression for a 
sphere? being 


1-—a=(6/2?) > (1/n?) exp (—n*x*Di/R?), (3) 


where R is the radius of the sphere. 

It was pointed out that because of the ex- 
tremely rapid convergence of the series given 
above, especially in the case of the slab, one could 


* Now at the University of Pennsylvania. 

1R. Ellickson and R. Ward, paper presented at the 
April, 1941 meeting of the American Chemical Society, to 
appear in J. Am. Chem. Soc. 

2 An expression similar to this was evidently first ob- 
tained by Fourier to represent the mean temperature of a 
sphere placed in a constant temperature bath. 
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obtain values of D, or rather D/L?, by plotting 
log (1—a) vs. t and measuring the slope of the 
curve at large values of ¢, when it approaches a 
straight line. This procedure is quite satisfactory 
in many cases but suffers from one limitation. 
The plot will approach a straight line when all 
terms in the series become negligible compared to 
the first term. This occurs when a reaches a value 
of about 0.3 for the case of a slab or about 0.7 
for the case of a sphere. This means that if one 
is dealing with diffusion in spherical particles, the 
diffusion process must be more than seventy per- 
cent complete before one can use the slope of the 
curve mentioned above to obtain a value of 
D/R*. Much of the data that might be used to 
check the equations does not give values of a 
greater than about 0.5, however, and this repre- 
sents a serious limitation on this method. 

It is the purpose of this paper to give a simple 
method for obtaining a value of D/L? or D/R* 
from values of a however small. In fact one value 
may be obtained from each experimentally de- 
termined value of a. If one writes Eq. (1) in 


TABLE I. Values of a calculated from Eqs. (1) and (3). 











SLAB SPHERE 
kt a a2/kt a [1 —(1 —a)!/3] 2 /t 
0.05 | 0.161 0.516 0.226 0.133 
0.1 | 0.227 0.516 0.310 0.136 
0.2 | 0.321 0.516 0.421 0.138 
0.4 | 0.454 0.516 0.560 0.143 
0.6 | 0.554 0.512 0.652 0.147 
0.8 | 0.636 0.505 0.721 0.150 
1.0 | 0.702 0.493 | 0.774 0.152 
1.5 | 0.819 0.447 0.864 0.157 
2.0 | 0.890 0.396 0.918 0.160 
3.0 | 0.960 0.307 0.970 0.158 
5.0 | 0.995 0.198 0.996 0.141 
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DETERMINATION OF 


the form, 


a=1—(8/m?) >) (1/n*) exp (—n?kt), (4) 
n odd 

where k=2°D/L?*, it is clear that a value of ki 
uniquely determines a value of a and vice versa. 
It is not possible to express @ as a function of kt 
in terms of any simple functions. It can be ex- 
pressed as an indefinite integral of a theta func- 
tion* but it is doubtful that this is of much use 
for calculations. However, it is very easy to 
calculate a table of values and plot one variable 
as a function of the other. The results of such a 
calculation from Eqs. (1) and (3) are given in 
Table I and the curves obtained are shown in 
Fig. 1. These curves represent the contribution 
of all terms in the series and therefore can be 
made as accurate as desired. Now if one takes an 
experimentally determined value of a, the corre- 
sponding value of kf can be read from the curve 
of Fig. 1. Knowing the value of ¢ to which this 
value of a, and therefore kt, corresponds one can 
calculate kt/t=2°D/L? and therefore determine 
D/L?. This procedure can be repeated for every 
experimentally determined value of a, and a 
value of D/L? is obtained for every experimen- 
tally determined point. The same procedure can 
be used for the spherical case except that the 
other curve is used, of course. 
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Fic. 1. a as a function of kt from Eqs. (1) and (3). 
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_* Whittaker and Watson, Modern Analysis (Cambridge 
University Press, 1940), fourth edition, p. 464. 
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Fic. 2. Diffusion of oxygen in copper. D is expressed in 
cm? sec.~! 


DIFFUSION OF OXYGEN IN COPPER 


The method described above can be applied 
to the data of Ransley‘ on the diffusion of oxygen 
in copper. These data represent the measurement 
of the amount of oxygen that has diffused from 
a strip of copper 0.09 cm thick as a function of 
the time. The quantity a is then the ratio of the 
mass of oxygen evolved at any time to the mass 
originally present in the strip. Some of these data 
were treated in our earlier paper. Figure 2 shows 
a plot of log (1—a) vs. t and it is observed to give 
a fairly good straight line at large values of ¢; 
the intercepts lie quite near the theoretical value 
of 0.8. However, the data taken at lower tem- 
peratures cannot be handled by that method and 
we shall treat it by the method described in this 
paper. 

In Table II are given values of a, ¢, kt, and 
kt/t; the last of these should be a constant and 
the table shows this to be fairly true. The letters 
(A) and (B) indicate two different samples of 
copper with different initial concentrations of 


4C. E. Ransley, J. Inst. Metals 65 (1939). (Advance 
copy.) 








oxygen, the values being 0.41 and 0.24 percent, 
respectively. The fact that sample B had the 
lower initial concentration may account for the 
fact that it seems to check better with the theory. 
This may indicate some dependence of D on 
concentration, a possibility that we have not 
considered here. The table shows that the data 
agree quite well with the theory. Figure 3 shows 
a plot of Eq. (4) with the average value of k as 
obtained above. It is probably true that the 
agreement is within experimental error. 

Ransley obtained rather good agreement with 
experiment using the parabolic law a?=ct; one 
might expect this since this law is derived from 
an assumption that is a fairly good approxima- 
tion to that upon which Fick’s law is based. The 
assumption, that the rate of diffusion through a 
thin strip is inversely proportional to the thick- 
ness of the strip, is an especially good approxima- 
tion for small thicknesses. It is an interesting fact 
that the series of Eq. (1) can be quite well repre- 
sented by the parabolic law in the early stages of 
the diffusion process. Table I shows that the 
parabolic law is equivalent to the exact solution 


TABLE II. Data on diffusion of oxygen in copper. 
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t (HR.) a@ kt kt/t 
600°C (A) 
62.5 0.077 0.0115 1.84 10-4 
110.0 0.097 0.0187 1.70 
205.0 0.131 0.0331 1.61 
700°C (A) 
2.0 0.0547 0.0051 2.55 X 107% 
4.5 0.0824 0.0132 2.93 
7.0 0.0960 0.0180 2.57 
8.6 0.1100 0.0236 2.76 
12.0 0.1164 0.0273 2.27 
800°C (A) 
0.5 0.0907 0.0156 3.12 X10 
1.0 0.1300 0.0328 3.28 
2.0 0.1713 0.0585 2.92 
3.5 0.233 0.104 2.97 
5.0 0.254 0.124 2.48 
8.0 0.308 0.182 2.28 
12.0 0.377 0.280 2.33 
800°C (B) 
0.25 0.075 0.011 4.40 x 10 
0.5 0.098 0.020 4.00 
2.0 0.202 0.078 3.90 
3.0 0.256 0.125 4.17 
4.0 0.298 0.170 4.25 
5.0 0.330 0.210 4.20 
7.0 0.383 0.286 4.10 
9.0 0.436 0.374 4.15 
12.0 0.500 0.490 4.10 
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Fic. 3. Diffusion of oxygen in copper. Circles are experi- 
mental points. Solid lines are theory from Eq. (4). 


within one percent for values of a up to about 
0.6; for values of a greater than this the ap- 
proximation becomes rapidly poorer. 

The advantage in obtaining a value of the 
diffusion coefficient as outlined in the present 
paper is that the value so obtained is directly 
comparable with those obtained in many other 
diffusion experiments, especially self-diffusion,° 
where there is little doubt but that Fick’s law 
holds. 

The temperature dependence of D is given by 
the equation, 


D=A exp (—Q/RT), (5) 


where Q is the activation energy, R is the gas 
constant and T the absolute temperature. It is 
evident that a plot of log D vs. 1/T should give a 
straight line with slope equal to Q/R. Such a plot 
is shown in Fig. 4. The activation energy deter- 
mined from the slope of the line is 50.5 kcal. per 
g-atom (2.19 ev per atom) as compared with the 
value of 46.0 kcal. per g-atom obtained by 
Ransley from a similar plot of the data treated 
according to his theory. 


5 F. R. Banks, Phys. Rev. 59, 376 (1941) and references 
there cited. 
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Fic. 4. Temperature dependence of diffusion coefficient 
of oxygen in copper. Activation energy from slope is 
50.5 kcal. per g-atom. 


APPLICATION TO SOLID PHASE REACTIONS 


The rate of many chemical reactions that take 
place in the solid phase is probably governed by 
the rate of diffusion of one of the reactauats. 
Assuming the particles to be spherical one can 
apply the theory developed here to such reac- 
tions. The data of Jander® on the rate of evolu- 
tion of carbon dioxide in the powder reaction 
between barium carbonate and silicon dioxide at 
830°C are given in Table III, together with 
values of kt and kt/t as determined from the curve 
of Fig. 1 for spherical particles. It is evident that 
the data agree with the theory remarkably well 
except at the very beginning of the reaction. The 
lack of agreement here is very probably due to 
the fact that there is a time lag before the entire 
mixture reaches a steady temperature. If one 
allows one minute for this time lag the agreement 
is excellent over the whole range. Jander esti- 
mated this time at two minutes but the correc- 
tion is necessarily quite approximate. The re- 
mainder of the data on this reaction taken at 
various temperatures from 800°C to 890°C agrees 
equally well with the theory. The values of kt/t, 
which are proportional to D/R?, are plotted on a 





*W. Jander, Zeits. f. anorg. allgem. Chemie 163, 1 
(1927). 
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logarithmic scale against 1/7 in Fig. 5 and it may 
be seen that the points lie quite close to the 
expected straight line. The value of the activa- 
tion energy (which does not depend on R) is 46.6 
kcal. per mole as compared with the value of 
56.0 kcal. per mole obtained by Jander. 

Data from the same author on the reaction 
between calcium carbonate and molybdic oxide 
at 493°C are given in Table IV. Again it is clear 
that the agreement with theory is excellent. 
Figure 6 shows a plot of log (kt/t) vs. 1/T for this 
reaction at various temperatures from 487°C to 
550°C. The value of the activation energy, 81.3 
kcal. per mole, obtained from the slope is some- 
what lower than the value of 86.3 kcal. per mole 
given by Jander. 

Using the same assumption that leads to the 
parabolic law for the linear case, that material 
diffusing through a layer diffuses at a rate in- 
versely proportional to the thickness of the 
layer, Jander derived an expression for the rate 
of reaction for a spherical particle, 


[1-(1-a)! P=ct, (6) 


where ¢ is a constant depending on the size of 
the particle and the diffusion constant. The last 
column in Table I shows that this is not nearly 
as good an approximation to the exact solution 
based on Fick’s law as is the parabolic law for 
the linear case. 

The fraction of completion of the diffusion 











TABLE III. Rate of evolution of CO2 in powder reaction. 
BaCO;+SiO:z 830°C 
t (MIN.) a kt (kt/t) X104 

3 0.021 0.0010 3.33 

4 0.031 0.0015 3.75 

5 0.040 0.0021 4.02 

6 0.047 0.0026 4.33 

7 0.054 0.0031 4.44 

8 0.061 0.0039 4.87 
10 0.071 0.0047 4.70 
15 0.090 0.0075 5.00 
20 0.100 0.0090 4.50 
25 0.108 0.0106 4.24 
30 0.120 0.0130 4.33 
40 0.134 0.0165 4.12 
50 0.146 0.0200 4.00 
60 0.163 0.025 4.16 
70 0.176 0.030 4.28 
80 0.187 0.033 4.13 
90 0.195 0.037 4.11 
100 0.205 0.041 4.05 
110 0.220 0.047 4.26 
120 0.227 0.050 4.16 
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barium carbonate and silicon dioxide. Activation energy is 
46.6 kcal. per mole. 


process at any time evidently depends very 
strongly on the size of the particles of the diffus- 
ing substance. The value of kt/t obtained for a 
reaction with different sizes of particles at the 
same temperature should indicate to what extent 
the theory is valid in this regard. The value so 
obtained should be inversely proportional to the 
square of the radius of the particles used. In 
Table V the value of D/R? from Jander’s data 
on the reaction between barium carbonate and 
silicon dioxide at 860°C for four different sizes of 
particles is given. Multiplying this value by the 
corresponding value of R? should give a constant 
value for D; the last column in the table shows 
the result to be quite satisfactory. 

It should be emphasized that the method of 
determining D outlined in this paper is applicable 
only if one is dealing with particles of uniform 
size. If this is not the case, the solution that is 
valid is not that of Eq. (1) or (3) but a sum of 
such series, one for each size of particle. In this 
case the method outlined here is quite useless. 
One can still use the earlier method of plotting 
log (1—«a) vs. t and the slope of this line at large 
t will be proportional to D/R®, where R is now 
the radius of the largest sized particles present. 
With a continuous distribution of particle size, 
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Fic. 6. Temperature dependence of reaction rate of 
calcium carbonate and molybdic oxide. Activation energy 
is 81.3 kcal. per mole. 





the problem would be very difficult to solve, 
unless perhaps the distribution could be repre- 


TABLE IV. 








CaCO3:+Mo0O3; 493°C 
a k 




















1 (MIN.) 1X10? (kt/t) X104 
2 0.019 0.10 5.00 
3 0.030 0.15 5.00 
4 0.042 0.22 5.50 
5 0.050 0.29 5.80 

10 0.078 0.57 5.70 
15 0.093 0.80 5.33 
20 0.103 1.00 5.40 
30 0.125 1.43 4.77 
40 0.152 2.20 5.50 
50 0.171 2.8 5.60 
60 0.185 3.3 5.50 
70 0.198 3.8 5.43 
80 0.213 4.4 5.51 
90 0.226 5.0 5.65 
100 0.234 5.3 5.30 
TABLE V. 
BaCO3+SiO2 860°C 
R (cm) D/R? D X10" (cm? sEc.”') 
0.0153 8.711078 2.04 
0.0086 2.06 X 1077 1.52 
0.0053 5.401077 1.52 
0.0036 1.16 10-* 1.50 
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sented by some very simple function. In any case, 
if there is more than one size of particle, the 
intercept will lie at a value lower than the value 
0.6 predicted by Eq. (3). 

In conclusion, one can say that the equations 
that kave been developed seem to have some 


4 


7 





experimental verification. The method outlined 
here greatly extends the type of data that can 
be compared with the theory. Further experi- 
mental data, especially that covering a diffusion 
process carried nearly to completion, are desirable 
as a further check. 
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Thermodynamic Properties of Binary Solid Solutions on the Basis of the Nearest 
Neighbor Approximation 


EpwIn N. LASSETTRE AND JOHN P. Howe* 
Laboratory of Chemistry, Ohio State University, Columbus, Ohio 


(Received August 28, 1941) 


An approximate method of computing the partition function of a binary solid solution is 
formulated. If the only constraints are on the mean energy and mean composition, it is shown 
that no metastable phase is predicted. The partition function is shown to be obtainable from 
the largest eigenvalue of a quadratic form and the condition for a phase transition is shown to 
be related to the degeneracy of the largest eigenvalue. The physical interpretation of the 
eigenfunction is shown to be related to the probability of surface configuration, while the 
square of the eigenfunction is related to the probability of a configuration on the interior of 
the crystal. Some simple examples are discussed which are related to the effect of coordination 


number on phase transitions. 


HE purpose of this paper is to develop a 
method of calculation of the approximate 
partition function for a binary solid solution. It 
is assumed that molecules of the two components 
are distributed at the lattice points of a given 
lattice and that the energy of a configuration is 
determined by the interaction of nearest neigh- 
bors. The properties of this model have been 
studied by many investigators with the use of 
various approximate methods.! Since the critical 
points predicted are of considerable interest, we 
feel that a further investigation will be useful. 
The statistical problem has been formulated 
by several people? but, to avoid confusion as to 


_—— address, Brown University, Providence, Rhode 
Sland, 

‘(a) W. L. Bragg and E. J. Williams, Proc. Roy. Soc. 
A145, 649 (1934); (b) H. A. Bethe, Proc. Roy. Soc. A150, 
552 (1935); (c) J. G. Kirkwood, J. Chem. Phys. 6, 130 
(1938); (d) for a review and bibliography see F. C. Nix 
and W. Shockley, Rev. Mod. Phys. 10 (1938); (e) R. H. 
Fowler and E. A. Guggenheim, Proc. Roy. Soc. A174, 189 
(1940); (f) H. A. Bethe and J. G. Kirkwood, J. Chem. 
Phys. 7, 578 (1939); (g) J. G. Kirkwood, J. Chem. Phys. 
8, 623 (1940). 

* The statistical problem for solutions has been pre- 
sented by E. A. Guggenheim, Proc. Roy. Soc. A135, 181 





standard states, a brief presentation is included. 
At the same time, certain general properties of 
the model are obtained. 


I. FORMULATION OF STATISTICAL PROBLEM 


For the purpose of numerical calculation, it is 
useful to formulate the statistical problem in 
terms of Gibbs’s grand ensemble. In terms of this 
scheme, the probability of finding 1, particles of 
type 1 and m2 particles of type 2 in a volume V, 
and state j is given by 


> Q+ 1x1 +N2x2— Ejnins 
P(j, 11, N2) =exp 
kT 





where £;,,n, is the energy of the jth state of a 
system containing m, particles of one kind and n» 
of the other. Q, x1, x2 are constants to be deter- 


(1932), and more recently the same author has discussed 
the application of the grand ensemble to the problem of 
regular solutions [Proc. Roy. Soc. A169, 134 (1938) ]. 
3’The principles underlying the grand ensemble are 
clearly presented in the book Principles of Statistical 
Mechanics by R. C. Tolman (Oxford, 1938), pp. 613-648. 















































mined by the conditions 


2, Pj, my m2) =1, (1) 
F ys mP(j, 11, M2) =i, (2) 
2, maP(j, m1, 2) = Ma. (3) 


Here, 7; and “2 are the mean numbers of par- 
ticles of each component. is interpreted as — PV 
where P is the external pressure and V is the 
volume of the system. x: and x2 are the chemical 
potentials of components 1 and 2, respectively. 
The work content of the system is 


A= —PV+i1x14+Nox2. 


In the present problem we consider a lattice 
consisting of M equivalent points and of such an 
extent that the total volume is Mv. The molecules 
of components 1 and 2 occupy these lattice points 
and a molecule of type 1 can replace a molecule 
of type 2 without any volume change or other 
distortion of the lattice. 

The energy of a state depends, in general, on 
the composition and is assumed for our purposes 
to have the form 


Ejinyny = Eoj+ 2%(2,E1+n2E2) +72V 2, 


where Eo; is the vibrational energy and depends 
on the size of the lattice but not on the arrange- 
ment of the particles, and E;, E2, and Vie are 
constants. The number of contacts between un- 
like molecules in a given configuration is 7, and 
x is the coordination number of the lattice under 
consideration. The microscopic state of the sys- 
tem is then given by specifying the vibrational 
quantum state j, the numbers of molecules 1; and 
m2, and the positions of the molecules in the 
lattice. A definite arrangement of the molecules 
in the lattice will be called a configuration. In 
specifying a configuration, like molecules are to 
be treated as indistinguishable. An exchange of 
two or more like molecules does not lead to a 
new configuration. 

As reference states the pure substances will 
be used so we consider the conditions (1), (2), 
(3) for the case 7;= M, fi2=0. Accordingly, the 
left side of (3) must vanish, and since each term 
in the sum is positive we must have P(j, m1, m2) =0 
unless m2=0. This can be approached by giving 
x2 a large negative value. We shall confine our 
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attention to the case for which this is exactly 
satisfied. The only configuration is, then, the one 
for which m.=0 and hence for which 1=0. We 
thus obtain 


—P\V+Mx°—3MxE; 
kT 


where [Ro ]“=2; exp (—Eo;/kT). When M is 
sufficiently large, Ro will approach to a constant 
independent of M and in what follows we sup- 
pose that M is so large that this is the case. P;° is 
the vapor pressure of pure component 1 and x,° 
is the chemical potential of the pure component. 
Using (4) we find 


xE,/2=x1°—Pyv+kT In Ro, (5) 





CRo]™ exp| |- 1, (4) 


where v is the volume per lattice point, v= V/M, 
and according to our assumption is independent 
of composition. In a similar way 


xE2/2=x2°—P2v+kT In Ro. (5’) 


We now introduce the following variables 








pact a 
a,=exp ; 
kT 
v(P2°—P)+x2—x2° 
a2=exp | (6) 
kT 
A\=a;/a2; Z2=exp {-— Vi2o/RT}, 
and Eqs. (1) and (2) become 
ag M= Le Ds v(M, Nn, a)d"2" (7) 
and 
Myx Jae 
—— — =. (8) 
ae Or 


v(M,n,i) is the number of configurations in 
which there are M lattice points, n molecules of 


type 1, and 7 contacts between unlike molecules. . 


Equations (5) and (5’) have been used to elimi- 
nate E,; and £2. 

The change in work content on mixing 7% 
molecules of component 1 with #2 molecules of 
component 2 is 


A =n kT log aytnokT log Qo. 


The condition for phase equilibrium in a 
heterogeneous system can be expressed in terms 
of a; and a». Consider, for convenience, a system 
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consisting of two phases, I and II, and let a;! 
and a;'! be the values of a in the two phases. In 
order that the two phases be in equilibrium in 
the absence of external fields we must have 
a} =a," and a2! =a,!!. These functions are thus 
similar to but not identical with activities. 

Certain useful properties of this model can 
be established without approximation and with- 
out complicated analysis. In the interests of 
clarity, these properties and the proofs are given 
at this point. 

(a) The function a for a component increases 
monotonically as the mole fraction of that com- 
ponent increases. 

We note that (1/Aa2)” decreases as \ in- 
creases, as is obvious from Eq. (7) when we 
remember that the number of configurations 
must be positive. Also, using Eqs. (7) and (8), 
we find that 


On, 
A— = ((1) wv) — ((11)w)?, 
Or 


where ((;")4,) and ((1))? are, respectively, the 
average of the square and the square of the 
average number of molecules of type 1. The 
derivative of 7%; with respect to X is, therefore, 
positive, and hence 7%, increases as \ increases. 
Since a,;=aeA, it follows that a; increases with 
increasing A, and hence also with increasing 7). 

(b) ¥(M, n, 1) =¥(M, M—n, i). 

Consider a configuration of molecules of 
type 1 on M lattice points with 7 contacts. If we 
replace each molecule of type 1 by a molecule of 
type 2 and vice versa we obtain a new config- 
uration of M—n molecules of type 1. The num- 
ber of contacts between. unlike molecules is 
unchanged. By means of this process, we obtain 
a one to one correspondence between configura- 
tions with m molecules of type 1 and those with 
M—n molecules of type 1, and it follows that 
V(M, n, i)= ¥(M, M—n, i). 

(c) ag 1(d) = Nag (1/A). 

From Eq. (7) 


as™(A) =A” +. i VM, 2, 1)(1/d) *—*2! 
= Y yn Di V(M, M—n, i)(1/d)*—"2". 


Using (b), the right side of the equation is just 
Mao (1/). 
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(d) #(1/A) = M—A(A). 

This can be proven directly from (8) by means 
of a straightforward differentiation. 

(e) a;(7) =a2(M— fi). 

This follows immediately from (b), (c), and (d), 
together with the fact that \a2=a1. 

As a corollary, we note that the change in 
work content on mixing is symmetrical around 
the point 7;=3M. This fact has been previously 
pointed out by Shockley.‘ It is to be noted from 
(d) that when A\=1, i= _M/2 and hence the mean 
composition is determined for this special value 
of X. It is also of interest that a is a monotone 
function of the mole fraction. This fact shows 
that in the present model we obtain no solution 
corresponding to a metastable phase unless some 
special constraint is imposed. It is of interest to 
see just how the phase separation enters in the 
subsequent treatment. 

In order to obtain an exact value of the func- 
tion as” it would be necessary to find the 
numbers ¥(M, n, 7), an exceedingly difficult task 
which has not been attempted here. Instead, we 
deal with the function a.~™” directly and find as 
a result that exact computations are possible in 
certain simple cases. 


II. DERIVATION OF THE DIFFERENCE 
EQUATIONS 


The partition function satisfies a set of linear 
difference equations with constant coefficients.*® 
This may be obtained in the following way. 

We consider a simple cubic lattice, although 
the method may be extended without difficulty 
to any type of lattice. The cubic lattice is as- 
sumed to consist of N parallel planes, each con- 
taining m lattice points. The lattice points are 
occupied by m molecules of type 1 and (mN—n) 
molecules of type 2. An arrangement of molecules 
in a given lattice plane will be called a distribu- 
tion. In the first plane we may have any one of 
2” distributions, and we denote by ¥;(N, 1, 7) 
the number of configurations for which the dis- 
tribution in the first plane is the /th distribution. 
If the first plane has the /th distribution, the 


‘W. Shockley, J. Chem. Phys. 6, 130 (1938); «bid. 6, 
523 (1938). 

5 This derivation follows the same general argument as 
that used by Fowler and Rushbrooke in treating a similar 
problem. Trans. Faraday Soc. 33, 1272 (1937). 
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second plane must have one of 2” distributions, 
say the kth. If 7; denotes the number of contacts 
between unlike molecules in the /th distribution 
and v,, the number of contacts of unlike mole- 
cules between the /th and kth (assumed ad- 
jacent), then 


VAN, n, i) =D Vi(N—-1,n—pi,t—ti—vn), (9) 


where p; is the number of molecules of type 1 in 
the distribution /. The sum extends over all 
possible distributions in the second plane. Multi- 
ply Eq. (9) by 2‘A" and sum over all ” and 7. We 
thus obtain after rearranging 


&(N) =2'r”: Do, 2d, (N—1), (10) 
where 
©(N)=don Di B1(N, n, 2)2°d". (11) 


Equations (10) are a set of 2” simultaneous linear 
difference equations with constant coefficients. 
When @,(1) is given, the solution is determined 
uniquely by these equations. 

The function ,(1) is the partition function for 
the /th distribution in a single lattice plane. This 
function plays the role of an initial condition, 
and we may consider it as being, to some extent, 
arbitrary. If, for example, the first layer of mole- 
cules is formed on the surface of an adsorbing 
crystal, $,(1) will depend on the nature of the 
crystal. With ®,(1) arbitrary, we have the prob- 
lem of finding asymptotic solutions of (10) which 
account for the thermodynamic properties of a 
large lattice. It is shown in the following two 
sections that the solutions are obtainable if we 
impose certain reasonable requirements on the 
asymptotic forms of the partition function. 


III. REDUCTION TO ALGEBRAIC FORM 


The grand ensemble has the property that the 
entropy of a homogeneous macroscopic phase 
under a given set of conditions is proportional to 
the quantity of the phase.* The same is true for 
the energy and hence for the work content. This 
latter condition is satisfied provided the quan- 
tities x; and x2 become independent of Nm as N 
and m become large. Thus, from Eq. (7), we see 
that 6(N)=2,#,(N) is the Nth power of a num- 
ber u(u=a2-”") as N becomes large, i.e., ®(V) ~u%. 
Moreover, the ratio ;(N)/(N) =¢: is the prob- 
ability of finding the /th distribution in the first 
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lattice plane, and as the number of lattice planes 
becomes large this number is expected to ap- 
proach a constant value independent of NV. These 
considerations suggest that we search for solu- 
tions which have the asymptotic form 


&(N)~om. (12) 
To find such solutions, substitute ¢,u* into (10) 
and divide by »*— thus obtaining 


udi=Zir?t Do, 2g; (13) 


Equation (13) is a set of 2” homogeneous linear 
equations in 2” unknowns ¢, and the condition 
of compatibility is 


det. | siuyPig’te — 161). | =(), (14) 
where 
(0, 1k 
bu= < 
\1,/=k 


and the eigenvalues are roots of Eq. (14). Since 
the secular equation (14) is of degree 2”, 2” roots 
exist. We designate the rth such root by uy, and 
the corresponding solution of (13) by @,. In the 
usual way, the following orthogonality condition 
is found to be obeyed 


> $1rb1,/Z"A"'=0 if Mr Us. (15) 


In what follows, we assume that the solutions 
have been normalized in such a way that 


D1 bir? /2trr'=1. (16) 
Equations (13) arise when the quadratic form 
S= or De Bikar TW}. (17) 


is made an extreme by varying zm; subject to the 
constraint 


Yi r2/siidrr=1. (18) 


This is obvious if we use Lagrange’s method of 
multipliers.6 The values of 7, which make (17) 
an extreme are solutions of (13), and the corre- 
sponding value of the quadratic form S is a root 
of (14). 


6 The general maximum and minimum properties o! 
quadratic forms is discussed in Methoden der Mathemati- 
schen Physik by R. Courant and P. Hilbert (Springer, 
1931), pp. 1-38. 
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IV. INTERPRETATION OF SOLUTIONS 


There remains the problem of selecting the 
eigenvalue of Eq. (13) which is of physical 
interest. Since ¢,/2.:¢: gives the probability of 
distribution /, it must always be positive. We 
now show that if yw: is the largest root of Eq. (14), 
and if the root is distinct, then the corresponding 
solution ¢ has the property ¢::/2:¢2 0 and is 
the only solution having this property. 

As previously stated, the extremes of the 
quadratic form S are the eigenvalues of Eq. (13). 
The largest eigenvalue must then be the maxi- 
mum value which S can attain. Consider a set of 
numbers z; which satisfy (18) and which have the 
property that 7,>0(/¥7z), 7;<0. The terms of S 
which are negative are 2”'i2,7;(147). If we vary 
the set by replacing 7; by —2; these terms be- 
come positive, and the magnitude of S increases. 
This variation obviously satisfies the constraint. 
Hence if the numbers 7, make S a maximum, 7; 
must have the same sign as 7,(1#7). In a similar 
fashion, it can be shown that 7, has the same 
sign for all values of /. The above argument shows 
that the best trial set is one in which the numbers 
are either all positive or all negative. Hence the 
solution $1, corresponding to the largest eigen- 
value, has either the property that ¢;,2 0 (all 2) 
or én <0 (all J). 

Now consider any other solution @,,. Since the 
largest eigenvalue is distinct, ¢), is orthogonal to 
¢n as in Eq. (15). This is only possible when ¢,, 
changes sign as / changes. 

The above argument demonstrates that if we 
substitute any arbitrary real numbers 7, in (17) 
which satisfy the constraint (18), we obtain a 
lower bound to the eigenvalue of physical inter- 
est, ie., S< yu. This result is useful for the purpose 
of making approximate numerical calculations. 
In the future, we designate $1 simply by ¢:. 

Using the fact that 6(N)~y* and the formula 
for mean composition previously given, we find 
that 

Nd Op 


——=f. 


pw Or 


This follows since, from (7), u is equal to a2”. 
By means of a straightforward calculation shown 
in the appendix, we find that 


N p pide /2'h?1= n. (19) 
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From this formula, the mean composition can 
be calculated directly from the solutions to 
Eqs. (13). 

It is interesting that Eq. (19) implies (but does 
not prove) that ¢//z'A”: gives the probability 
of configuration / in a lattice plane far in the 
interior of the lattice. A proof of this is included 
in the appendix. 

We have thus far discussed the case in which 
the largest root is distinct. When the largest root 
is degenerate, we have several linearly inde- 
pendent solutions, any linear combination of 
which is also a solution, and hence the surface 
probabilities are indeterminate. In case the roots 
cross, certain properties of the system will, in 
general, change discontinuously, since the largest 
root is always the one of physical interest. It 
seems natural then to regard the temperature at 
which the largest root becomes degenerate as a 
transition temperature. 

The foregoing argument shows that a sufficient 
condition for a transition is that the largest root 
becomes degenerate. In case the transition is a 
phase separation and mm is finite, this condition is 
also necessary. This follows, since the condition 
for phase separation is that at two mean com- 
positions %, and wi», the values of a@ for com- 
ponent 1 must be equal, and similarly for 
component 2; therefore, A\.=, and pa=ps. Since 
(NX/u)(dp/dA) = 7, it follows that 

Ou) Op 
—| #—, 
Or, 


and the first derivative of »% with respect to A is 
discontinuous. However, uw is the root of an 
algebraic equation with analytic coefficients, and 
hence du/0X is continuous unless the largest root 
is degenerate. This latter statement is valid only 
when m is finite. 

We can easily see how the composition is af- 
fected from Eq. (19). By taking linear combina- 
tions of the solutions corresponding to the de- 
generate roots, we obtain from (19) a range of 
compositions rather than a single composition. 
In other words, the activity is independent of 
composition over a finite range of compositions. 

From the results of Section I the phase 
boundary must be symmetrical around mole 
fraction 4 and hence i; =mN—/A, and \.=A=1. 
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This means that a degeneracy which leads to 
phase separation occurs when X is unity.’ 

On the basis of numerical calculation for 
special cases, the authors feel that the largest 
root becomes degenerate (at temperatures above 
absolute zero) only in the limiting case as m 
. approaches infinity. It is true, however, that as 
the largest and next largest roots approach, the 
system tends to simulate one in which a transi- 
tion is taking place. 


V. SOLUTION OF SPECIAL CASES 


The general Eq. (13) has not been solved in 
closed form, but certain special cases can be 
solved analytically, and these form the basis for 
solution by perturbation theory. Simple cases 
which can be treated numerically are also of 
interest, since they can be compared directly 
with the approximate calculations of other 
investigators. 

Consider first the case of an ideal solution, i.e., 
one for which z=1. The right-hand side of Eq. 
(13) is then independent of / and hence if 1.¥0 


o1=AX”,, (12a) 


where A is independent of /. A can be determined 
from the normalizing condition (16) the result 
being A =(1+A)”/?. If we substitute in (13) the 
value of u is obtained 


w=(1+A)”. (13a) 
This is the only non-vanishing eigenvalue and 


7 An exception to this may occur if the phase boundary 
curve (in the temperature-composition diagram) has a 
minimum at mole fraction 3. This possibility is not rig- 
orously excluded by the results of Section I, but it seems 
rather remote and is ignored in what follows. 
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hence is the largest. It follows immediately from 
(7), and the fact that (N) ~y, that ag=1/1+4. 
From Eq. (19) we find that the mole fraction of 
component 1 is X¥,;=A/1+A and hence a,;= Xj, 
a2=X >. These conditions are not those ordinarily 
obtained for perfect solutions but are in agree- 
ment with the considerations recently advanced 
by Bury.* When the vapor pressures of the pure 
components are equal, these conditions reduce 
to Raoult’s law, i.e., the activity of a given com- 
ponent is equal to the mole fraction of that 
component. 

Consider next the case in which m=1. This 
corresponds to a linear lattice.® 

The equations are 


Loi = 1 +292, 
20 
(u/A)b2=261+- G2. - 
The roots of the secular equation are 
2u=1+A+((1—A)?+4Az2*)!. (21) 


These roots become equal only when z=0 and 
\=1. When A <1, the largest root is 1 and i=0. 
When A>1, w=A and *=N. When A=1, w= 1 
and from Eq. (19) we see that the mole fraction 
may vary between 0 and 1, the activity being 
independent of composition and equal to unity. 
This illustrates the behavior of the equations at 
a transition point. 

Using Eqs. (20) and (21), a can be computed 
as a function of composition; the results are 
shown in Fig. 1 for the case z=}. 

A valuable contribution to the theory has been 
made by Kirkwood,! who showed how to com- 
pute the coefficients in the expansion of the work 
content in powers of 1/kT. In a recent paper, 
phase separation was discussed, and the condi- 
tion for occurrence!® in terms of our notation 
was found to be 


— Vie/kT. = (1—4/x)!—1. (22) 


When 7 <T,, a separation into phases can take 
place. We note first that the coordination number 
x must be equal to or greater than 4 in order that 
T. be real. It should be mentioned in this con- 
nection that other theories of order and disorder 


8 Bury, Trans. Faraday Soc. 36, 795 (1940). 


9 A rigorous discussion, equivalent to this, has been given 
by Ising, Zeits. f. Physik 31, 253 (1925). 
10 J. G. Kirkwood, J. Phys. Chem. 43, 97 (1939). 
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also lead to results which are similar but not 
identical with this. For example, with the Bragg 
and Williams approximations x21, with the 
Bethe approximation x23, while in the recent 
quasichemical method of Guggenheim and 
Fowler!© x23. The investigation of a simple 
lattice of coordination number four has, there- 
fore, a certain interest. 

A simple lattice having coordination number 4 
is that for which m=3, and the lattice planes 
contain three points at the corners of an equi- 
lateral triangle. In this case, there are 8 distri- 
butions of which two sets of three have equal 
probabilities by symmetry. We designate the 


Oo 
distinct configurations by the scheme o o (1) 
x Oo x 

o o (2) x x (3) x x (4). The simultaneous equa- 
tions are then 


Ubi = $1+3262+32"o3+2% 4, 
uds/2’A = 261+ (14+22")b2+ (22+2°)b3+2°%o, 
uds/2°d? = 271+ (22 +2°) h2+(1+22") os +264, 

uds/d3 = 2361 +32°h2+ 3263+ ou. 


We have solved these equations numerically for 
z=4 (i.e., Vie/RT=+1.0986). The results are 
shown in Fig. 2. It is evident from the figure that 
although a is not constant it changes very slowly 
with mole fraction around mole fraction 3. When 
\=1, the secular determinant factors into two 
quadratics. The roots of these. quadratics are 
never equal for positive values of z except in the 
trivial case z=0. On this basis it seems that some 
property of the lattice other than coordination 
number must be involved in order to assure a 
phase separation. It is undoubtedly true that 
extension of Kirkwood’s expansion to higher 
powers of Vi2/kT will introduce other properties 
of the lattice than coordination number and thus 
distinguish between the three-filament lattice 
and an extended two-dimensional lattice with a 
coordination number of four. In fact, Kirkwood 
and Bethe'“) have recently shown that, for 
lattices which are made up of two interpenetrat- 
ing lattices, other properties than coordination 
number enter in the coefficient of (Vi2/RkT)*". 





" As a check on the correctness of Eq. (13) we have 
expanded A in powers of Vi2/kT for the case of a square 
plane lattice with coordination number four, and with 
A=1. This was done by expanding z, u, ¢ in powers of 
Vi2/kT, the coefficients being determined by substituting 
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The method developed here is complicated by 
the fact that algebraic equations of high degree 
must be solved. It has been our experience that 
the process of solution, especially near a critical 
point, involves elaborate analysis. This analysis 
appears to be inevitable, since the mathematical 
behavior of the thermodynamic functions is not 
simple in this neighborhood. 

This complication is offset, in part, by the fact 
that the eigenfunctions have a direct physical 
significance and only the largest eigenvalue is 
necessary. It is pleasing also that the criterion for 
a critical point can be stated very simply in terms 
of degeneracy of the largest eigenvalue. Finally, 
the existence of a variational principle is of some 
assistance in correlating and refining approxi- 
mate calculations. Progress toward obtaining 
numerical solutions has been made and will be 
reported in a later paper. 


APPENDIX 


I. Proof of equation for the mean number of 
particles, Eq. (19) 


We have previously shown that 
(1a) 
(2a) 


B= Doi Die B*hide, 
> or siyr=l, 


If we differentiate each of these equations with 
respect to A, we obtain, taking account of the 





1.0 U T I 4, | 1 
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in (13). The expansion coefficients thus obtained agree 
with those of Kirkwood and Bethe (reference 10) if we set 
s equal to zero in the latter and substitute the proper values 
of the coordination number and ¢. The fifth moment is 
also in agreement with that recently reported by T. S. 
Chang, J. Chem. Phys. 9, 169-174 (1941). 
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set of simultaneous equations satisfied by ¢i, 
Ou 09): 
—=2n >> (o— /rnsi), (3a) 
Or k Or 
O91 —— 
21> (6 /rva) => (pigs? /d"'z"t) (4a) 
l Or l 
and hence, eliminating 


dg 
pe (+ rs) 
l Or 


Ou n 
—-—=D (pid? /d?8") =—. (Sa) 
porn t N 


we get 


This is Eq. (19) which we wished to prove. 


II. Probability of a configuration in a lattice plane 
far within the interior of a lattice 


Consider a lattice plane in the interior of the 
lattice. We suppose that there are /—1 planes 
to the left and N—1 planes to the right, and that 
the particles of the plane are arranged in accord- 
ance with the /th distribution. The partition 
function for the assembly is, using our previous 
notation and maintaining the /th distribution: 


Q:=6,(M)@,(N) /ds". (6a) 


We divide by z'A”: because, by definition, each 
of the functions ®,(17) and ®,(N) contains this 
product as a factor, and hence the numerator of 
(6a) contains the square of z2''\?:. The probability 
of distribution / on the interior lattice plane is 
then 


Q:/>1 Qi. 


If the lattice plane is sufficiently far in the in- 
terior of the lattice, we may replace ®,(M) by 
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om”, &(N) by dm and we obtain 
12/dP ist (7a) 


for the probability, taking account of the normal- 
ization of ¢:. 


III. Asymptotic form of solution 


The results of Sections III and IV may be 
developed in an alternative way by using an 
exact solution of Eq. (10). To obtain this solution 
we expand #,(JV) in a series of solutions of (13). 
The expansion coefficients may be obtained by 
substituting in (13) and using the orthogonality 
relations. The final result is 


©(N) =D, Curr. 


C, is a constant (i.e., independent of NV) which is 
obtainable from 


1 (1) 1, 


br § ZtiyPa 





r 


If u; is the largest eigenvalue and is distinct 
&(N) = Cidnmr®(1+O(exp {| —N log (u:/u,)}) J. 


It is immediately clear that when JN is sufficiently 
large &)(N)~Cidnui®. Equation (12) is thus 
justified. In case the largest eigenvalue is de- 
generate, ®,(V) is obviously a linear combination 
of two or more solutions. In an actual case for 
which m is of order N? we may encounter a dis- 
tribution of eigenvalues for which log (u1/u,) is 
of order 1/N for a range of values r. We have 
not been able to examine this possibility in a 
general way. However, in the limiting cases for 
which exact solutions are obtainable (i.e., at 
z=0, 1, etc.) the eigenvalues are well behaved. 
This fact leads us to believe that log (u:/y,) is 
not of order 1/N except in the neighborhood of 
a real transition point. In a later paper this 
question is investigated in some detail for one 
case. 
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The Thermal Decomposition of Methyl-n-Butyl Ether! 


SIDNEY J. MAGRAM AND H. Austin TAYLOR 
New York University, New York, New York 


(Received July 3, 1941) 


The thermal decomposition of methyl-n-butyl ether proceeds in a manner analogous to that of 
the other ethers. There is no evidence to support the suggestion that it might decompose into 
methyl alcohol and ethylene by a ring closure through a hydrogen bond. In the first stages of 
the decomposition about 60 percent of the ether decomposes into CH, and C;H;CHO and 
about 40 percent into HCHO, C,H, and C;:H;. The decomposition of the aldehydes occurs at 
rates greater than that of the ether, probably induced by radicals from the ether. The butyralde- 
hyde yields mainly CH4, C2H, and CO. First-order constants fall off as the reaction proceeds and 
also with diminishing pressure. The energy of activation at about 200 mm pressure is 54.0 
kcal. and at 100 mm 56.6 kcal. For the reaction inhibited by nitric oxide the energy of activa- 


tion is about 60 kcal. 


INTRODUCTION 


N a recent paper Rice and Teller? have ex- 
amined the question whether product mole- 
cules of a decomposition are formed in a single 
elementary act or through a chain mechanism, 
from the point of view of the principle of least 
motion. According to this principle the activation 
energy and total motion of atoms required to 
bring the atoms from the old to the new position 
should be as small as possible. On such a basis it 
appeared that butyl] derivatives offered the possi- 
bility of forming ring structures with subsequent 
direct rearrangement into product molecules. 
Thus for methyl butyl ether, decomposition 

might be expected as follows: 

CH, 

, 

HsC—O CH: 
: | —CH;0H+2C>H,. 

H CH, 


NZ 
CH 


Besides presenting a difference in reaction prod- 
ucts from those found from other ethers, the 
above reaction, were it to proceed, should occur 
at lower temperatures than those required for 
the chain decomposition, a simplification not to 
be overlooked in view of the present meager ex- 
perimental data on unimolecular reactions. The 
occurrence of the reaction depends on the O—H 





' Abstract from a thesis submitted by S. J. M. in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy at New York University, February 1941. 

*F.O. Rice and E. Teller, J. Chem. Phys. 6, 489 (1938). 
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“hydrogen bond”’ that must form. Configura- 
tionally the methyl hydrogen comes close to the 
oxygen atom, especially at higher temperatures 
when bending vibrations are excited. It appeared 
plausible, therefore, that at least in part the de- 
composition might occur by this path and de- 
terminations of the methyl alcohol produced 
reveal the quantitative extent of this reaction. 


EXPERIMENTAL 


Commercial methy! butyl ether with a boiling 
range of 69-72°C was fractionally distilled yield- 
ing a distillate between 70.1-—70.3°C at 750 mm 
pressure. This was again distilled in contact with 
sodium and at low pressure from room tempera- 
ture to —80°C using solid COs, several times, 
being finally stored over sodium at — 80°C. 

A liter Pyrex reaction vessel, enclosed in an 
electrically heated furnace, was attached to a 
constant volume manometer, ether storage bulbs, 
vacuum pumps, two Toepler pumps and a gas 
analysis apparatus. Since the vapor pressure of 
methyl butyl ether at room temperature is only 
about 12 cm it was necessary to heat the con- 
necting tubing from the reaction vessel to 
about 50°C to prevent condensation. As a lu- 
bricant in the stopcocks which were heated, a 
plasticizer, di-o-xenylmonopheny! phosphate was 
found satisfactory. 

In making a run, ether was allowed to vaporize 
into a mixing chamber surrounded by hot water 
to the desired pressure. In experiments with 
nitric oxide the latter was introduced first into 
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TABLE I. 
CONTRAC- O2 Con- 
Temp.°C | P/Pi | VOL. (cc) CO2 TION SUMED 
475 2.0 6.7 13.6 17.2 24.2 
482 1.7 2.0 4.2 5.3 7.5 

















this chamber. The gas mixture was permitted to 
expand rapidly into the evacuated reaction 
vessel. No attempt was made to equalize the 
pressures. The change of pressure with time was 
followed for the desired period and the products 
were withdrawn by one of the Toepler pumps and 
their volume noted. The gases were then passed 
back and forth through a trap surrounded by an 
ice-toluene or liquid-nitrogen bath between the 
two Toepler pumps until the pressure was ob- 
served to drop to a constant value. At dry ice 
temperature undecomposed ether and _inter- 
mediate aldehydes were condensed. All other 
gases, even higher possible hydrocarbons, would 
pass through at the pressures employed. 


NORMAL DECOMPOSITION 


Analytical results 


The material condensed in the gas trap showed 
reducing properties on Tollens reagent. A frac- 
tion of a cc gave a silver precipitate from a solu- 
tion of AgNO; in NH,OH and NaOH. 

One cc of the condensed material heated with 
5 cc of a ten percent NaOH solution gave a 
yellowish brown resinous material characteristic 
of the lower aliphatic aldehydes. The resorcinol- 
sulfuric acid test showed formaldehyde present. 
A still more definite identification was obtained 
by heating a few tenths of a cc with a trace of 
chromotropic acid in strong sulfuric acid. A deep 
violet color is shown by formaldehyde. Acetal- 
dehyde, propionaldehyde and butyraldehyde do 
not give this test. 

To establish the identity of the other aldehyde, 
which qualitative tests showed was present in 
larger amounts than formaldehyde, a derivative 
was made by heating a few cc of condensed ma- 
terial with a saturated alcoholic solution of 2.4 
dinitrophenylhydrazine. This when warmed with 
a drop of concentrated HCI gave a yellow pre- 
cipitate upon addition of water. The melting 


°F. Feigl, Qualitative Analysis by Spot Tests. 
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point of the compound after one crystallization 
from alcohol was 120-121°C. It showed no de- 
pression when mixed with the hydrazone of 
butyraldehyde. Since 147°C is the nearest melt- 
ing point of a hydrazone, that from acetaldehyde, 
it is certain that the aldehyde here found is 
butyraldehyde. 

Although many tests were made for methyl 
alcohol, no evidence for its presence, even in 
small amounts, could be found. 

Analysis of the uncondensed gases showed the 
presence of methane, unsaturates, carbon mon- 
oxide, hydrogen and a higher hydrocarbon. The 
percentage of unsaturates was considerably larger 
than one would expect merely from analogy with 
the products of the lower ethers and the em- 
pirical formula of the saturated hydrocarbon 
mixture was considerably lower in carbon than 
was expected. Thus by analogy with the over-all 
decomposition of the lower ethers the following 
products might be expected: 


HCHO+ C,H 
7 —-H2+CO+CiHi 
CH;0CH.,CH2CH2CHs 


N\ 
CH,+C;:H;CHO 
—CH,+C3;H;+CoO. 


Only small amounts of unsaturates produced 
from the hydrocarbons and a total hydrocarbon 
content of average formula C2Hg is to be ex- 
pected. Results show that the amounts of un- 
saturate are higher early in the reaction than 
later, that is, they constitute about 30 percent 









































TABLE II. 
CONTRAC- O2 Con- 
Temp. °C | P/Pi | Vor. Gas | COz TION SUMED 
477 3.32 13.1 26.8 30.9 44.6 
489 1.97 20.8 41.7 47.2 68.1 
533 1.85 20.0 28.6 44.0 62.6 
535 2.45 9.4 19.1 20.8 30.5 
TABLE III. 
Temp. °C VoL. Gas CO2 ConT. Oz 
477 0.49 1.00 1.15 1.67 
489 50 1.00 1.14 1.64 
533 2 1.00 1.14 1.62 
535 49 1.00 1.09 1.60 
C2He 50 1.00 1.25 1.75 
C2H, 50 1.00 1.00 1.50 
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DECOMPOSITION OF METHYL-BUTYL ETHER 


TABLE IV. 











% CONDENSED % IN 


P/P: Temp. aT —80°C KOH 


2) 
° 


REMARKS 





480 
510 
520 


496 
509 
509 
515 
496 
497 
506 
500 
464 
520 


519 
535 


2.97 542 
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of the uncondensed gas so that the over-all formula 
of the hydrocarbon mixture is as low as C;.3H,4.7. 

Such a relatively high amount of unsaturate 
might, at first sight, suggest that the decomposi- 
tion into methyl alcohol and ethylene was oc- 
curring, at least in part. The complete absence 
of methyl alcohol, the lowness of the percentage 
of unsaturates (even if the methyl alcohol were 
completely decomposed into CO and Hz, the 
ethylene percentage should still be 40), the 
amounts and complexity of the other products 
and finally the similarity of the temperature 
range in which methyl butyl ether decomposes to 
that of the lower ethers all argue against the 
ultimate molecule split. 

To identify higher hydrocarbons they were 
condensed in liquid nitrogen and pumped free of 
H2, CO and CHy. Complete combustion after 
removal of unsaturates gave the data of Table I 
for the products of two runs withdrawn at the 
temperatures and ratios of pressure to initial 
pressure indicated. The ratios of the last four 
columns are 0.49 : 1.00 : 1.26: 1.78 and 0.48 : 
1.00 : 1.26 : 1.79 at 475° and 482°C, respectively. 
The ratio for C.H¢ is 0.50: 1.00 : 1.25 : 1.75 
showing that, within the experimental error for 
the small volumes used, the gas is essentially 


ethane. It should be noted that CsHs and CyHyo 





. are fairly stable at these temperatures or de- 


compose very slowly. 

The identity of the unsaturate was then ob- 
tained from a complete combustion of the un- 
saturate and ethane mixture. Several such results 
are given in Table II. In the form of ratios these 
data appear as in Table III. The ratio of volume 
of gas to CO, indicates that the unsaturate must 
be only C,H,y. The values of contraction to 
oxygen consumed are consistent with mixtures 
of C,H, and C;He¢, though one cannot rely ac- 
curately on these values to give the ratio of C.H, 
to C.H¢ since, during the pumping, more C,H, 
than C,H. would be lost. 

Table IV records a number of analyses at 
various temperatures and various stages of de- 
composition. The small absorption in KOH may 
be due to CO». 

From these analyses and the intermediates 
identified it is clear that at least two simultaneous 
processes are occurring. One would appear to be 
a reaction yielding formaldehyde and hence 
C,H as a residue. This residue must give C,H, 
and C.,H,. A radical mechanism can account 
for this. 


R+CH;0C,H,—-RH+CH20C,Hg, 
CH.OC,H,—-HCHO+C,Hsg, 
C,H »—-C2H,+ CoH s. 
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‘ 4 all, — ) 
7) 20 se #0 
Time 
RuN Temp. °C Pi (Mm) TIME UNIT (MIN.) 
1 543 56 2.10 
2 510 59 9.80 
3 511 84 7.72 
4 484 99 28.0 
5 560 125 0.91 
6 504 111 8.25 
7 494 115 13.0 


Fic. 1. P/P; against time. The unit of time for each run 
is the time when AP=P;/2. 


The over-all reaction would be 
CH;,0C,H»—>HCHO+C.Hi+ CoHe. (1) 


Similarly the formation of butyraldehyde by a 
radical mechanism might be represented : 


R+C,H,OCH;—-RH+CH;CH.2CH2CHOCHs;, 
CH;CH»CH2,CHOCH; 
—CH;CH.CH,CHO+CHs3. 


This over-all reaction would be 


CH;0C,H»—CH,+C;3;H;CHO. (2) 


Since at the beginning and well on in the reac- 
tion CH, is in excess of CsHe, it appears that 
C,H, must also be produced from butyraldehyde. 
This, however, should be accompanied by more 
CO than appears in the analyses. For example, 
at 515° where P/Pi=1.8, 30.1 percent C2H,, 16.3 
of CO, 9.3 of He and 16.9 of C2H¢ are found. 
Assuming that H, and CO come principally from 
formaldehyde in reaction (1), then 16.3—9.3=7.0 
percent CO must have come from reaction (2). 
However, since in reaction (1) C2.H, and CzH,g are 
produced in equivalent amounts, then 30.1— 
16.9=13.2 percent C2.H, must have been pro- 
duced from reaction (2). This excess of C2H,4 over 
CO suggests that butyraldehyde decomposes 
first into C.H, and a residue containing oxygen. 
Specific tests for acetaldehyde gave negative 
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results, as also did the indigo test for compounds 
containing the CH;CO group. Furthermore, Rice‘ 
found no evidence of ketene or acetaldehyde as 
intermediates in butyraldehyde decomposition. 
It would appear likely that butyraldehyde de- 
composes as does acetaldehyde 


R+C;H;,CHO—RH+C;3H 7CO. 














This latter radical would be less stable than 
acetyl and split CO readily leaving the propyl 
radical which is known’ to break down into C2H, 
and CHs, the latter finally yielding CH,. The 
excess of C.H, over CO can be partially ac- 
counted for on the basis of an observed hetero- 
geneity of about 10 percent of the whole reaction. 
Davoud and Hinshelwood found a heterogeneous 
decomposition of diethyl ether yielding ethylene 
and water as opposed to the homogeneous reac- 
tion yielding acetaldehyde and ethane. 

A balance of reaction products suggests, there- 
fore, that from 58 to 63 percent of the ether de- 
composes giving butyraldehyde, the remainder 
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i i 1 
‘ 
Tine ” . 
RUN Temp. °C Pi (wm) TIME UNIT (MIN.) 
1 475 176 25.2 
2 489 243 10.7 
3 513 285 3.72 
4 527 276 1.75 
5 536 207 1.43 
6 541 240 1.14 





Fic. 2. Log (P:—P;) against time. The unit of time for 
each run is the time when AP = P;/2. 






4F.O. Rice, private communication. 

5 Rice, Walters, and Ruoff, J. Chem. Phys. 8, 262 (1940). 

6 Davoud and Hinshelwood, Proc. Roy. Soc. A171, 39 
(1939). 
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giving formaldehyde. The decomposition of these 
intermediates is fast compared with the initial 
ether decomposition suggesting decomposition 
induced by radicals from the ether. 

No accurate balance of products can be made, 
since some material was formed just outside the 
furnace in the connecting tube to the reaction 
vessel. It is not known whether this originates 
from the ether or from its decomposition prod- 
ucts. Furthermore, the analyses of the products 
show quite definitely that the percentages of 
hydrogen and ethylene decrease while that of 
ethane increases as the reaction proceeds; hydro- 
gen is reacting with ethylene. 

It is interesting to note that the products found 
and accounted for by a radical mechanism re- 
quire reaction of the radical with a hydrogen 
atom on carbon adjacent to the oxygen, sub- 
stantiating previous observations that C—H 
bonds in @ positions (i.e., 8 toa C=O group) are 
more reactive. If hydrogen from the 1, 2, or 3 
carbon atoms in the butyl group had been re- 
moved, percentages of hydrogen from 25 to 40 
and unsaturates higher than ethylene would be 
expected. 


Pressure change data 


Although interpretation of the kinetics of ether 
decompositions from pressure-time curves is al- 
most impossible it is nevertheless advisable, for 
comparison with other ethers, to note the simi- 
larities found in methyl buty] ether. A glance at 
Fig. 1 shows immediately that the slower reac- 
tions subsequent to the initial ether decomposi- 
tion make the final pressure doubtful. To com- 
pare runs at different temperatures on the same 
plot, the time unit for each run is taken as the 
time required for the pressure to increase by fifty 
percent. Attempts to estimate where the reaction 
appeared to be going if decomposition of the 
intermediates did not occur led to ratios of 
calculated final to initial pressure varying con- 
siderably with temperature and pressure, for 
example, from 2.8 at 120 mm to 3.1 at 60 mm. 
On the other hand, ratios of actual final to 
initial pressure varied up to 4.0, though these 
values cannot be other than approximate as the 
pressure is still slowly rising. 

A plot of log (P:— P;) against time for different 
temperatures in Fig. 2 shows that the relative 
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ETHER 
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Fic. 3. Inhibition by NO. 





rates of initial and secondary decompositions up 
to fifty percent decomposition of the ether do not 
change much over the temperature range studied. 
The rate falls off slightly faster at lower tem- 
peratures. The times for the pressure to increase 
25, 50, and 100 percent of the initial pressure 
correspond roughly to 12.5, 25, and 50 percent of 
the ether decomposition. For a first-order reac- 
tion these times should be in the ratio of 0.19 : 
0.45 : 1.2. The average values found for a large 
number of runs at pressures varying from 100 to 
to 200 mm and temperatures from 460° to 560°C 
were 0.18 : 0.40 : 1.00 indicating that the reac- 
tion is sensibly first order.? However, calculating 
the velocity constant k; from the initial slopes of 
the pressure-time curves, a plot of log k; against 
log P; has an initial slope of 0.5, indicating an 
over-all rate of 3 order, in the range from 25 to 
about 100 mm pressure. Above 200 mm the 
reaction is approximately first order. 

The energy of activation calculated from, ki, 
4; and #4 at pressures from 10 to 100 mm averaged 
56.6 kcal. At pressures above about 200 mm the 
value was 54.0 kcal. In each case a small increase 
in the activation energy 55.7 to 57.3 and 52.9 to 
55.5 in going from k; to ti2, that is as the reaction 
proceeds further, was observed. This may be 
attributed, as in similar cases,® to a chain inhibi- 
tion by some intermediate product such as 
ethylene. 

It should be remarked that the above data 
could be reproduced only by adopting a uniform 
procedure preliminary to a run. Prolonged pump- 
ing for several hours always produced a de- 


7 The actual data are contained in the original thesis of 
S. Magram in New York University Library. 

8 Fletcher and Rollefson, J. Am. Chem. Soc. 58, 2131 
(1936); F. O. Rice and O. L. Polly, J. Chem. Phys. 6, 273 
(1938). 

















760 Ss. j. 
TABLE V. Inhibition of methyl butyl ether decomposition 
: by NO. 
THE N 
Temp. °K — aad Rinhib Rinhib/Fnorm 
758 90 0.03 0.017 0.90 
773 95 13.50 0.007 0.14 
773 105 6.40 0.006 0.13 
773 97 3.20 0.006 0.13 
773 103 1.40 0.013 0.28 
773 99 0.13 0.019 0.40 
776 128 6.70 0.007 0.10 
786 81 0.12 0.038 0.43 
787 90 11.10 0.011 0.12 
790 225 3.50 | 0.033 0.20 
793 283 2.30 0.042 0.17 
793 280 4.30 0.040 0.17 
793 275 0.15 0.180 0.74 
798 158 2.30 0.046 0.20 
803 120 1.70 0.053 0.23 
803 255 8.00 0.061 0.18 
807 102 2.30 0.040 0.19 
808 253 0.80 0.127 0.30 
809 106 3.20 0.042 0.18 
809 240 0.30 0.191 0.46 
809 245 8.00 0.080 0.20 
812 102 0.25 ° 0.086 0.31 
813 79 1.50 0.078 0.31 




















composition showing a high initial rate which 
was not reproducible. Pumping for from five to 
ten minutes only, after a run, enabled subsequent 
ones to be repeated. To investigate further this 
obvious surface effect, the reaction vessel was 
partly filled with glass tubes with fire polished 
ends whereby the surface to volume ratio was 
increased five times. When reproducible runs 
were obtained with this increased surface, the 
rates were an average of 1.45 times faster than 
the normal rates. This indicates that about 10 
percent of the normal reaction is heterogeneous 
in nature. 


Inhibition by NO 


The inhibited decomposition of methyl butyl 
ether was similar to that found by Hinshelwood 
and Staveley® in the lower ethers. Figure 3 shows 
two typical runs together with the uninhibited 
rates for comparison. The effect can be seen to 
be transitory from the acceleration of the rates. 
Whether NO is used up by the time the inflection 
in the curve is reached is not known. A test for 
cyanide made on the products collected from 
about eight runs was positive. This can be 


( _— and Staveley, Proc. Roy. Soc. A159, 192 
1 x 
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accounted for by decomposition of formaldoxime 
which is isomeric with CH;NO. Table V contains 
the data for the inhibition. These data show that 
the effect of from 3 to 13 mm NO on the ratio of 
inhibited and normal rates is constant. This is 
probably due to a balancing of a chain inhibition 
and a catalytic effect by NO following a mecha- 
nism suggested by Taylor and Burton for acetal- 
dehyde.’® This is further borne out by the ob- 
served decreased over-all! energy of activation 
with increased pressure as predicted by this 
mechanism. 

A plot of log kinni» against 1/7 for a few runs 
with initial ether pressures from 225 to 280 mm 
gives an energy of activation of about 60 kcal. 
The mean chain length shows a small decrease 
with increased temperature. This would manifest 
itself as an increased temperature coefficient of 
the inhibited reaction. 

Table VI presents comparative data for a 
number of ethers. The constancy of the energies 
of activation of both normal and inhibited de- 
compositions shows that the same mechanism is 
probably responsible in all cases. 

There is no evidence at these elevated tempera- 
tures of a hydrogen bonded ring formation. Since 
this might occur at room temperature, however, 
a preliminary photolysis” was attempted mainly 
to see if ethylene is produced in considerable 
amounts. The situation would be analogous to 
that found by Thompson and Dainton" in the 
photolysis of the alkyl nitrites. The ether vapor 
was exposed in a quartz reaction vessel to the 
radiation, mainly 2537A, from a low pressure 
quartz mercury arc, spiral in form and _ sur- 
rounding the reaction vessel. A pressure increase 











TABLE VI. 

ETHER MEAN CHAIN LENGTH Einhib Enormal 
CH;0CH; 17.0 at 540°C 62.0 58.5 
CH;0C2H; 8.4 at 540°C 62.0 54.5 
CH;0C;3H;7 3.1 at 545°C 49.01 —_ 
CH;0C,H, 5.0 at 536°C 60.0 - §4.0 








10 Taylor and Burton, J. Chem. Phys. 7, 414 (1939). 
11 This low value is probably due to catalysis by iodides. 
See Ure and Young, J. Phys. Chem. 37, 1183 (1933). 
2 The authors wish to thank Dr. A. Gordon of this 
laboratory for this work. 
( 13 — and Dainton, Trans. Faraday Soc. 33, 1546 
1937). 
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PYROLYSIS OF FORMALDOXIME 


indicated that decomposition occurred. The 
fixed gases in the products were withdrawn by 
means of a Toepler pump through a trap in 
liquid nitrogen. A relatively considerable volume 
was so withdrawn. Withdrawal of a second frac- 
tion was attempted at the temperature of melt- 
ing sec-butyl chloride which should contain the 
ethylene. The amount of this was less than 0.1 
cc. This result was checked in a second run and it 
was decided that the split into methyl alcohol 


761 


and ethylene was not occurring and further work 
was abandoned. 

It appears then that the principle of least 
motion, in predicting a direct rearrangement into 
product molecules of large molecules with ex- 
posed negative groups, must be restricted to such 
cases where, for example, an electronic resonance 
exists between initial and intermediate states, 
adding to the strength of the new bond or bonds 
in the intermediate. 
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Formaldoxime decomposes in a first-order reaction to give hydrogen cyanide and water at a 
rate k=109-5e—39000/RT sec, The products react together at a slower rate to give principally 
carbon monoxide and ammonia although nitrogen, hydrogen, ethane and traces of unsaturate 
are also produced. The reaction is probably heterogeneous on a clean glass surface. It is acceler- 
ated by traces of air and by nitric oxide. These facts offer an adequate explanation of the fate 
of methyl radicals in nitric oxide inhibited reactions by the isomerization of the CH;NO 


complex into formaldoxime. 


INCE the time of the original observation? 

that small amounts of nitric oxide were able 
to inhibit almost completely the decompositions 
of numerous organic materials, the mechanism 
of the action has been the subject of much 
experimentation but of more speculation. It is 
quite generally admitted that the action in- 
volves a complex formation between a radical 
and nitric oxide, for example, with a methyl 
radical, CH3;NO. It is not so generally admitted 
that the complex is stable. In many pyrolyses 
inhibited by nitric oxide, tests for cyanides in 
the reaction products have been obtained. In 
others the presence of ammonia has been demon- 
strated. In still other cases a white solid has been 
isolated and it has been suggested*® that this 
solid was formaldoxime. 





' Abstract from a thesis submitted in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy, 
New York University, June, 1941. 
ass and Hinshelwood, Proc. Roy. Soc. A154, 335 


*Linnett and Thompson, Trans. Faraday Soc. 33, 874 
(1937); Davis, Jahn, and Burton, J. Am. Chem. Soc. 60, 
10 (1938); Taylor and Jahn, J. Chem. Phys. 7, 470 (1939). 


On the other hand however, reactions have 
been followed in presence of nitric oxide and no 
positive tests for cyanides or ammonia were 
found‘ wherefore decomposition of CH;NO into 
methane, carbon monoxide, nitrogen and water 
was proposed. Further it has been suggested 
that the stability of the complex relative to the 
reactants forming it is small. Echols and Pease 
suggested that an equilibrium is set up between 
the complex, nitric oxide and radicals in butane 
pyrolysis even at 520°C. 

Although other interpretations may be offered 
for some of the observations that have led to 
the above suggestions, the variety of these 
suggestions indicates a lack of knowledge of the 
stability of possible isomers of the complex. 
For the complex CH;NO, the most obvious 
isomers are formamide and formaldoxime. Of the 
latter little is known at all of its decomposition 
products and according to Sherman® formamide 


4 Staveley, Proc. Roy. Soc. Al62, 557 (1937). 


5 Echols and Pease, J. Am. Chem. Soc. 59, 766 (1937). 
6 Sherman, Am. Chem. J. 20, 226 (1904). 


























H. A. TAYLOR 
TABLE I. 

T =412° T =385° T =376° T =350° 

TIME Pi =103.5mMM Pi=68.5MM Pi=81.0MM Pi =129.3 MM 

MIN. AP (mM) AP (mm). AP (mm) AP (mM) 
1.0 6.4 0.2 0.5 0.3 
) 10.5 0.6 0.9 _ 
2.0 13.8 1.4 15 0.8 
2.5 16.7 1.6 1.9 — 
3.0 19.5 2.2 2.4 1.0 
3.5 22.0 2.9 2.9 — 
4.0 24.9 3.6 3.8 1.9 
4.5 26.8 4.0 _ _— 
5.0 29.0 — 4.7 2.8 
5.5 31.0 5.4 — — 
6.0 33.0 5.9 $5 3.8 
7.0 34.5 6.9 6.5 4.8 
8.0 35.3 re —_— — 
9.0 39.2 — 6.5 
10.0 41.9 8.9 —- 
15.0 S12 10.8 9.8 








gives ammonia and carbon monoxide on distilla- 
tion and when heated above 195°C yields also 
hydrogen cyanide and water. Whether isomeriza- 
tion from one to the other is possible is not 
known. 

The first attempt to follow the mechanism of 
nitric oxide inhibition directly was made by 
Thompson and Meissner.’ They followed the 
change in the absorption spectrum during re- 
action of a mixture of nitric oxide and dimethyl] 
ether. It was found that the nitric oxide bands 
progressively weakened as ether decompositions 
progressed and that ammonia bands appeared. 
The disappearance of the nitric oxide bands 
paralleled the resumption of normal ether decom- 
position. It was suggested that formaldoxime 
might be the intermediate. Using diethyl ether 
in a similar experiment new bands appear which 
could be attributed to methylamine and the 
following reaction was proposed : 


CoH;+ NO—CH;NH2+ CO. 


Alternatively an unstable nitroso compound may 
have been present since, when the reaction 
products were frozen out a blue solution resulted. 
The production of carbon monoxide might how- 
ever result from formamide, the corresponding 
ammonia reacting with radicals to give methyl- 
amine. The evidence therefore is not unequivocal 
and it was hoped that a study of the decomposi- 
tion of formaldoxime might throw light on the 
whole subject. 


7 Thompson and Meissner, Nature 139, 1018 (1937). 
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EXPERIMENTAL 


Formaldoxime was prepared by the method of 
Scholl. Hydroxylamine hydrochloride was dis- 
solved in a minimum amount of water, com- 
pletely neutralized with sodium carbonate and 
cooled below 5°C in an ice bath. Formaldehyde 
solution is added with constant stirring, the low 
temperature being maintained. A heavy white 
precipitate soon separates. The mixture is allowed 
to stand in an ice box for a few hours to ensure 
complete reaction and filtered on a Buchner 
funnel. The solid is washed with water at about 
50°C until free of chlorides and dried in air. The 
product may be purified by sublimation below 
100°C without decomposition. The white solid is 
a polymer of formaldoxime. 

To determine the molecular state of the vapor 
from this solid, the molecular weight was deter- 
mined by observing the pressure at which a 
known weight of the polymer as vapor occupied 
a predetermined volume at the chosen tempera- 
ture. The results of three such determinations at 
about 100°C were 46.9, 45.4 and 46.0. The 
molecular weight of CH» : NOH is 45.0, hence it 
may be concluded that the vapor is completely 
depolymerized at 100°C. 

In preliminary bomb tube experiments the 
presence of hydrogen cyanide, ammonia, and 
carbon monoxide was established among the 
products of decomposition. The bomb tubes after 
heating at about 400°C for varying lengths of 
time showed on cooling heavy brownish deposits 
which were attributed to polymerization of hy- 
drogen cyanide and drops of liquid which ,per- 
sisted on standing. The monomeric formaldoxime 
is a liquid under ordinary conditions but poly- 
merizes on standing. The liquid therefore was 











TABLE II. 

Temp. °C k X104 sEc.7! Temp. °C k X104 sec. 
413 10.8 385 3.0 
412 1.5 2 | 381 2.4 
412 10.1 378 2.4 
412 10.2 378 2.3 
408 10.0 376 au 
408 8.8 375 2.1 
405 15 374 1.9 
401 6.6 370 1.4 
398 5.1 365 $2 
386 2.9 350 1.2 











8 Scholl, Ber. 24, 573 (1891). 
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taken to be water. A small amount of unabsorb- 
able gas was always observed. It was decided to 
study the pyrolysis using a static method by 
following the pressure change during decom- 
position making subsequent analyses of the 
products. 

A 500-cc Pyrex reaction vessel in an electrically 
heated furnace was connected to a constant 
volume manometer, and through a mercury cut- 
off to the formaldoxime reservoirs, diffusion 
pump system and to a Toepler pump. To pre- 
vent condensation, all the connecting tubing was 
heated to about 100°C. As a lubricant on the 
heated stopcocks a plasticizer, di-o-xenylmono- 
phenylphosphate was used, the reaction system 
during a run being closed by mercury. The solid 
polymer of formaldoxime was vaporized into a 
liter storage bulb kept in boiling water. The 
reactant was drawn from this bulb as required. 

In the early experiments before the mercury 
cut-off was introduced and before the plasticizer 
was used as lubricant very erratic results were 
obtained. These were traced in part to leaks 
due to the heated stopcocks. The accelerating 
effect of traces of air may be significant if acting 
as chain starters. On the other hand a clean 
glass surface had an accelerating effect. The rate 
of reaction decreased steadily to a reproducible 
limit as successive reactions are carried out. 
Alternatively this condition could be obtained 
by decomposing acetone for a prolonged period 
in the reaction vessel prior to a run. The ad- 
mission of air to such a prepared surface com- 
pletely spoiled the reproducibility, subsequent 
rates being again high. The reaction is un- 
doubtedly in part heterogeneous on a clean glass 
surface though whether the poisoning of the 
surface by this method was complete was not 
tested by poisoning in like manner a more 
extended surface. 

At the conclusion of a run the mercury in 
the cut-off was lowered and the products with- 
drawn by means of the Toepler pump for 
analysis. 


RESULTS 


Experiments were made over the temperature 
range from 350° to 415°C. The ratio of final to 
initial pressure varied from 1.8 to 1.9. Analysis 
indicates that the principal primary products are 
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TABLE III. 

Temp. °C TIME Ve/Vm 
415 Zero 10.0 
415 2 days 3.1 
415 3 days y Be 
435 20 min. 8.0 
450 10 min. 7.7 
450 12 hr. 2.2 
450 3 days 1.3 





hydrogen cyanide and water and that a 2 to 1 
ratio for final to initial pressure should result. 
The most likely explanation of the lower final 
ratio is that hydrogen cyanide is a reactive sub- 
stance and especially in presence of water is 
known to polymerize. Analytical data presented 
later show that complicating secondary reactions 
of hydrogen cyanide and water become sig- 
nificant later in the over-all reaction. 

At the lowest temperatures studied the re- 
action shows a short induction period, the cause 
of which could not be found. At 370°C and above, 
however, the induction period is too short to be 
noticed and plots of the logarithm of the formal- 
doxime pressure against time are initially linear. 
The reaction interval during which linearity of 
this plot persists increases with temperature. At 
370°C the reaction deviates from first order after 
about one-fifth of the time for total decomposi- 
tion, while above 400°C the reaction is first order 
to almost the half-life. A possible explanation of 
the falling off in the rate and its variation with 
temperature again lies in the polymerization of 
hydrogen cyanide formed. At higher tempera- 
tures a relatively higher concentration of mono- 
mer will be present and the normal pressure 
increase will persist to a greater extent of decom- 
position of the formaldoxime. Alternatively at 
lower temperatures the higher proportion of 
polymer causes a less rapid total pressure in- 
crease and the reaction appears to fall off. In 
Table I are given some illustrative data of the 
observed pressure increase. 

Assuming the reaction to be the simple two- 
mole product for one of reactant, the pressure of 
undecomposed formaldoxime at any time may 
be taken as the initial pressure less the observed 
pressure increase. Velocity constants were calcu- 
lated from the slope of the best straight line 
through data of the logarithms of the pressure 
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against time. A summary of these constants is 
given in Table II. 

These data give a good straight line when 
plotted as log k against 1/T and the slope of the 
line corresponds to an energy of activation of 
39 kcals. 

The preliminary analytical work having shown 
the presence of hydrogen cyanide, ammonia and 


. water among the products of decomposition it 


was realized that these would condense or be 
dissolved in passage through the Toepler pump 
system and would not be collectible entirely for 
analysis. To estimate the approximate amount of 
such absorbable gases the products withdrawn by 
the Toepler were washed with dilute acid and 
alkali and their volume V,, measured. Knowing 
the volume of the reaction system the volume 
V. of gas, theoretically to be withdrawn, can be 
calculated from the pressure and temperature in 
the reaction vessel. Table III gives several 
values of the ratio of V./Vm at various: tempera- 
tures for the gases withdrawn at the times 
indicated. These times are measured from the 
time that reaction appears to be complete as 
judged by the pressure change. 

The ratio V./Vm is seen to decrease with time, 
that is, the proportion of gas not removed by 
acid and alkaline washing increases with time. 
If the only products of formaldoxime decom- 
position were hydrogen cyanide and water, 
V./Vm would be infinite. If however these can 
react further to give carbon monoxide and 
ammonia, V./V» for complete reaction would be 
unity. The observed trend of the ratio indicates 
that hydrogen cyanide and water may react 
further. This is borne out by some rough analyses 
of the washed gases which show that when the 
ratio V./Vm is high the percentage of carbon 
monoxide is only about 40. As V./ Vm decreases 
the percentage rises to over 60. The reaction is 
by no means simple since about 20 percent of 
the gas is nitrogen and small amounts are present 
of hydrogen, unsaturates and a hydrocarbon, the 
average ratio of contraction to carbon dioxide 
produced on combustion of which indicated it to 
be ethane. That the increase in production of 
carbon monoxide is probably not due to a di- 
rect decomposition of formaldoxime into carbon 
monoxide and ammonia simultaneously with 
decomposition into hydrogen cyanide and water 


is inferred from separate experiments. Equi- 
molecular mixtures of hydrogen cyanide and 
water were sealed in bomb tubes and heated for 
varying periods of time in a furnace at 450°C. 
The tubes were opened under water. The water 
on analysis showed the presence of ammonia as 
well as cyanide. The water only partly filled the 
tubes showing that insoluble gases had been 
produced. These on analysis showed an average 
of 61 percent carbon monoxide, 20 percent 
nitrogen, 15 percent hydrogen, 4 percent ethane 
and traces of unsaturate. These values are in 
close agreement with those found when the 
products from formaldoxime are allowed to re- 
main in the furnace for some time after pressure 
change had ceased. Finally to confirm that the 
nitrogen and other products are produced from 
hydrogen cyanide and water and not from carbon 
monoxide and ammonia, bomb tubes of the 
latter mixture were heated as previously for 
several days. On cooling a slight white deposit 
formed on the sides of the tube. On opening under 
water the deposit dissolved. The solution showed 
the presence of ammonia but no cyanide. The 
insoluble gas was almost pure carbon monoxide. 


DISCUSSION 


The pyrolysis of formaldoxime is seen to 
involve a primary dissociation into hydrogen 
cyanide and water with an accompanying pres- 
sure change. These products react further, but 
more slowly, to produce chiefly carbon monoxide 
and ammonia by a reaction almost devoid of 
pressure change. The velocity constants pre- 
viously calculated thus represent the rates of the 
primary dissociation and the activation energy 
calculated from them is for the reaction: 


CH,NOH>HCN +H.0. 
The rate of this reaction may thus be expressed: 


k = 109-5¢e—39000/RT sec.—. 


The low value of the frequency factor, indicating 
that the process of activation is accompanied by 
an appreciable decrease in entropy, is not un- 
expected in view of the restriction of motion 
that must exist in the activated formaldoxime 
when about to dissociate into hydrogen cyanide 
and water, over that in the normal molecule. 
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The value of the activation energy 39 kcal. is 
close to that found by Steacie® for the N—O 
bond in alkyl nitrites namely 37.7 kcal. This too 
indicates that it is the N—O bond in the ac- 
tivated complex which breaks. 

The subsequent reactions of hydrogen cyanide 
and water must remain for the present matters 
for conjecture. It is known!® that hydrogen 
cyanide forms several polymers. One of these, 
diiminoethylene, HN =C =C=NH might readily 
react with water to produce carbon monoxide 
and ammonia. Another, diamino maleic dinitrile 


H2N. NH: 


Mh 
C=C 


* \ 
NC 
might conceivably be a source of hydrocarbons, 
while the complete pyrolysis of higher polymers 
has been reported" to yield nitrogen and hy- 
drogen. 

The presence of ethane in the final products 
suggests a decomposition of formaldoxime into 
nitric oxide and a methyl radical. Many tests 
were made for nitric oxide and all proved nega- 
tive. Furthermore when nitric oxid was pur- 
posely added in very small quantities to formal- 
doxime, a marked increase in the rate of pressure 
change was observed. This increase was so 
marked that in one attempt to introduce nitric 
oxide when formaldoxime was already partially 

® Steacie and Shaw, J. Chem. Phys. 3, 344 (1935). 


10 Bedel, Bull. Soc. Chim. 35, 339 (1924). 
"St. Claire Deville, J. iiber die Forts. Chem. 307 (1863). 
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decomposed, a violent explosion ensued, shatter- 
ing the apparatus. The work was not pursued 
and no attempt was made to compare the 
products from formaldoxime in presence and 
absence of nitric oxide. The reaction, however, 
appears to be similar to that in presence of small 
amounts of air where also a marked acceleration 
of the rate of pressure increase is observed. If 
nitric oxide were produced during the course of 
the reaction an auto-acceleration would be ex- 
pected. The induction period at lower tempera- 
tures is the only observation of any such effect. 
Even if this was the correct explanation of the 
presence of ethane and of the induction period, 
the amount of formaldoxime decomposing in this 
manner must be small and the effect of the nitric 
oxide must be that of a chain inducer. The whole 
possibility seems doubtful in view of the forma- 
tion of ethane from hydrogen cyanide and 
water. 

The problem of the fate of methyl radicals in 
presence of nitric oxide can thus adequately be 
explained on the basis of the formation and 
subsequent decomposition of formaldoxime. The 
alternative production of formamide is not so 
satisfactory in accounting for many of the 
observations in nitric oxide inhibited reactions. 
Thus formamide is liquid whereas a white solid 
is usually observed. It decomposes readily to give 
ammonia and carbon monoxide at lower tem- 
peratures without any trace of hydrogen cyanide 
while at higher temperatures only small traces of 
cyanide can be found. Positive tests for cyanide 
are much more frequent than for ammonia in 
nitric oxide inhibitions. 
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An apparatus for the measurement of double refraction 
of flow is described. The magnitude of birefringence, ex- 
pressed as the rotation of the plane of polarized light, may 
be measured to 0.05° under the best conditions with the 
assembly. The vortical angle and double refraction of flow 
of dioxane solutions of a series of fractionated polymers of 
a-methyl methacrylate have been measured at 25° over a 
concentration range of from one gram per liter to 50 grams 
per liter, depending upon the polymer under investigation. 
The specific birefringence, referred to the viscosity of the 
solvent, is a constant characteristic of each polymer in 
sufficiently dilute dioxane solutions. With the possible 
exception of behavior in such dilute solutions, the vortical 
angle is dependent upon concentration. Values of the 
rotary diffusion constant are calculated for each polymer 


1. INTRODUCTION 


gee the many high molecular weight 
substances which have been investigated 
recently with the new techniques of colloid 
chemistry, a large number are found which 
possess elongated structural units. In dilute solu- 
tions the behavior of these macromolecules is 
described by the molecular kinetic theory. There 
are several procedures which may be followed in 
studying these dilute solutions to give informa- 
tion about the length and length to width ratio 
of the solute molecules. Three of them, sedi- 
mentation analyses,’ dispersion of dielectric con- 
stant,‘ and viscosity,> have been the subject of 
much theoretical and experimental development. 
Less well-grounded from a theoretical viewpoint, 
and at the present time more difficult for experi- 
mentation, is the measurement of the double 
refraction exhibited by such solutions when they 


1 More complete details of this work are to be found in 
the thesis of Armand J. de Rosset submitted to the faculty 
of the University of Wisconsin in partial fulfillment of the 
requirements for the Ph.D. degree in June, 1939. 

2 Present Address: Universal Oil Products Company, 
Riverside, Illinois. 

3 Svedberg and Pedersen, The Ultracentrifuge (Clarendon 
Press, Oxford, 1940). 

4P. Debye, Polar Molecules (Chemical Catalog Com- 
Page) New York, 1929); Perrin, J. de phys. et rad. 5, 497 

1 , 

5 Staudinger, Die hochmolekularen organischen Naturstoffe 
(Berlin, 1932); Mark, Physical Chemistry of High Polymeric 
Systems (Interscience Publishers, New York, 1940). 
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from measurements of the vortical angle taken in dilute 
solutions. The lengths of the polymers are calculated from 
the rotary diffusion constants. They vary from 890A to 
1880A. The shape and structure of the polymers are dis- 
cussed from the point of view of their electrostatic and 
steric properties as evidenced by a molecular model. 
Tentative values of molecular weight and of length to 
width ratio are assigned. The viscosity of dilute dioxane 
solutions of the polymers has been measured at 25°. The 
specific viscosity increment at infinite dilution is estimated 
by extrapolation. This quantity is a linear function of the 
square of the length to width ratio, in partial confirmation 
of several theories of the viscosity of macromolecular 
solutions. 


are subjected to a mechanical shearing stress. 
In this research the attempt was made to ob- 
serve the flow double refraction in dilute dioxane 
solutions of methyl methacrylate polymers and 
to estimate the size and shape of such polymers 
from these data. 


2. THEORY 


Double refraction studies are usually made in 
the Kundt cell. It consists of two concentric 
cylinders, one of which may be rotated at a 
uniform speed. The solution under investigation 
is placed in the annular space between the 
cylinders. A beam of plane polarized light is 
passed through this cylindrical shell of liquid in 
the direction parallel to the axis of rotation. When 
the solution is sheared, it exhibits birefringence, 
which may be measured according to the classical 
method as the apparent rotation of the plane of 
the polarized light, by using a quarter-wave 
plate compensator. 

In determining the birefringence of non- 
dichroic solutions, both the amount and direction 
of double refraction must be considered. If the 
plane of the incident polarized light is rotated 
through 360°, the magnitude of birefringence will 
pass through four maxima at intervals of 90°. 
One pair of maxima will make an angle equal to 
or less than 45° with the direction of flow. This 





tl 


Ss’ 


fc 


de 


ve 


al 
ex 





DOUBLE 


angle is known as the vortical angle and will be 
designated by the symbol y. The magnitude of 
the maxima is the amount of double refraction 
of the solution, and will be designated by the 
symbol B. From the theory of birefringence it 
follows that 

B=2Z/180h. (1) 


In this equation Z is the observed rotation of the 
plane of the polarized light in degrees, and \ and 
h are the wave-length of light used and the depth 
of medium traversed, respectively, measured in 
the same units. 

In this report use is made of a quantity S, 
defined as S=B/Gyp. In this expression G is the 
velocity gradient between the cylinders, and 1 
is the viscosity of the solvent. If the solvent 
alone shows a slight birefringence Bo a similar 
expression holds, 

So=Bo/Gno. 


Then the specific birefringence of the solute is 
defined as 
D=(S—So)/c, (2) 


where c is the weight concentration of solute in 
grams per liter. 

Practically all modern theories of double re- 
fraction of flow of solutions presuppose that the 
elongated solute molecules are oriented by the 
velocity gradient through the solution. The re- 
sulting birefringence may be attributed to the 
intrinsic optical anisotropy of the solute par- 
ticles; or it may arise from their shape alone, as 
was pointed out by Wiener.* The former type 
of double refraction is “‘intrinsic,’’ the latter is 
“form” double refraction. The amount of form 
double refraction depends on the refractive index 


of the solvent and, from Wiener’s theory, should _ 


vanish when this is equal to the refractive index 
of the solute particles. Sadron,’ however, deduces 
a somewhat different situation. 

The first step in developing an orientation 
theory of double refraction is to determine the 
angular distribution of the particle axes under 
the opposing influences of the velocity gradient 
tending to orient the particles, and of Brownian 
motion tending to disarray them. This has been 
done by Boeder,®’ who solves graphically the 

° Wiener, Physik. Zeits. 5, 332 (1904). 


”Sadron, J. de phys. et rad. 8, 481 (1937). 
* Boeder, Zeits. f. Physik 75, 258 (1932). 
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FLOW 


two-dimensional equation of forced diffusion, 
60°p/0¢°?+G9A(p sin? ¢)/d¢=0. (3) 


Here, @ is the rotary diffusion constant, p is 
the direction density of particle axes, and ¢ is 
the angle between a particle axis and direction of 
flow. This equation assumes that the particles 
are rigid, homogeneous, and that they are long 
enough to approximate liquid lines. It also pre- 
supposes that the particles do not interact to 
prevent independent rotary diffusion and that 
there is no other orienting influence but that 
of the velocity gradient. Boeder’s solution is 
given in Table I, where B and y are given as 
functions of the ratio G/@. The quantity B has 
been normalized to 27 as the maximum possible 
amount of birefringence. From the data in the 
third column, it is possible to determine immedi- 
ately the rotary diffusion constant of any system 
of macromolecules which satisfies Eq. (3) from 
an observation of the vortical angle at any given 
gradient. 

Kuhn,° using much the same hydrodynamic 
treatment which he employs in his theory of 
viscosity, relates the rotary diffusion constant to 
the length of the particle by the equation 


49 = Ono = 8k T/rs’, ( 4) 


in which @ is the rotary diffusion constant re- 
duced to unit viscosity of solvent mo, & is the 
Boltzmann constant, T is the absolute tempera- 
ture, and s is the length of the particle, assumed 
to be much larger than the width. Kuhn also 
relates the specific birefringence of the solute to 
the rotary diffusion constant by an equation 
which may be reduced to the simple form 


D6) =constant. (5) 


The constant is valid for a homologous series of 
macromolecules. One assumes that the total 
anisotropy of the macromolecules is proportional 


TABLE I. Boeder’s solution of the distribution problem. 








v / B 


45.00° 2.000 
41.45° 2.400 
38.25° 2.770 
32.22° 3.020 














® Kuhn, Zeits. f. physik. Chemie A161, 1 (1932). 
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to their mass. This is equivalent to assuming 
that the maximum potential birefringence of the 
solution (complete orientation of macromole- 
cules) is proportional to the weight concentration. 
Sadron’ has recently considered the effect of 
solvent, neglected by Boeder and by Kuhn. His 
general equation for double refraction flow is 


1 Béo No 
N Gno Ny? +2 
3 a*—b? 


= ae +A’ (eac1—es02) |, (6) 








where V=number of molecules per unit volume; 
a,b=major and minor axes of molecule, con- 
sidered as an ellipsoid of revolution ; 01, ¢2=opti- 
cal polarizabilities of the molecule along these 
aXeS ; a, €& = geometrical constants of the ellipsoid ; 
no= refractive index of the solvent ; A’ = (mo? —1)/ 
(no?+2). The term involving (o1—¢2) represents 
intrinsic double refraction, while that involving 
(€a¢1—€y02) represents form double refraction. 
The denominator of the left-hand side should be 
multiplied by a factor of 3 to correct for an 
algebraic error copied from an earlier paper by 
Langevin. We find, in attempting to reproduce 
Sadron’s integrations, that apparently an addi- 
tional factor of 4 should be introduced into the 
denominator as well. 

Equation (6) predicts that a solution of opti- 
cally isotropic macromolecules will always show 
double refraction of flow, whatever the solvent. 
That is, there will exist no solvent of critical 
refractive index in which the birefringence will 
vanish. This conclusion is inconsistent with that 
of Wiener who states that form double refraction 
(which is all that need be considered when 
dealing with intrinsically isotropic particles) will 
vanish when the solute and solvent have the 
same refractive index. Wiener’s theory is valid, 
however, only for particles much smaller than 
the wave-length of light. Many macromolecules 
do not satisfy this condition. 

Equation (6) may be reduced to a form com- 
parable with Eq. (5) by assuming as before that 
the anisotropy of the macromolecule is propor- 
tional to its mass. One obtains for a long aniso- 
tropic particle 


D « (1/60)[1—A’/A’. J, (7) 
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where A’, refers to the refractive index function 
for the critical solvent in which double refraction 
of flow disappears. 


3. APPARATUS AND METHODS 


For the viscosity measurements described in 
this report the Washburn-Williams modification 
of the Ostwald viscometer is used. The efflux 
time for water is 190 seconds at 25.00°C using a 
25 cc volume of liquid. The equation for kine- 
matic viscosity in stokes k, is 


k, =4.87 X10-5t—0.056/t, 


where ¢ is the time of efflux in seconds. Density 
measurements, required for the calculation of 
absolute viscosities, are made in a single bulb 
Ostwald pycnometer of volume nearly 8 cc. 

The apparatus used for measuring the double 
refraction of flow of polymeric methyl metha- 
crylate solutions consists essentially of a Kundt 
cell and a precision polarimeter mounted hori- 
zontally on an optical bench. 

The outer cylinder or stator of the Kundt cell 
is brass with an inner radius of 5.293 cm and a 
length of 10.71 cm. It has a permanently at- 
tached base plate and a cover plate which may 
be removed to permit insertion of the inner 
cylinder, or rotor. The base plate of the stator 
is provided with a recess which holds a bronze 
bearing and the rear end of the rotor shaft. 
The cover plate may be firmly clamped to the 
stator by a large threaded bronze ring. It is 
drilled through the center and equipped with a 
collar through which passes the front end of the 
rotor shaft. This collar houses the packing and 
packing nut, a Torrington No. T-201 ball race, 
and a large set screw which holds the ball race 
firmly against a shoulder on the rotor shaft. 
The packing used in the experiments on dioxane 
solutions was woven cotton wick material which 
had been soaked in a mixture of glycerol and 
300-mesh graphite. The cover plate is also 
equipped with an inlet tube and filling bulb, and 
a gauge to observe the level of the solution in 
the cell. 

Recessed slits, 3 mm long, are cut through both 
cover and base plates near their periphery to 
correspond to the gap between the cylinders. 
Brass window holders are held in the recesses by 
three screws. Windows cut from strain-free 
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optical glass are cemented into the holders with 
glycerol-litharge cement. When the window 
holders are assembled to the cell, the screws are 
tightened evenly to preclude any strains which 
would make the windows double refracting. The 
windows are aligned by a brass pin in the cover 
plate which drops into a hole to correspond in 
the rim of the stator. 

The stator is wrapped with 3” copper tubing 
through which water at 25°+0.03° from a large 
reservoir is driven by a circulatory pump. 

The rotor is constructed of phosphor bronze 
tubing, the brass end plates and steel shaft being 
soldered on. It has an outside diameter of 5.266 
cm and a length of 10.69 cm. The gap between 
the cylinders is 0.0268 cm, which is small enough 
to permit high velocity gradients to be attained 
before turbulent flow sets in. Clearance of the 
rotor at the base and cover plates of the stator 
may be adjusted by a set screw in the base plate 
which regulates the height of the bronze bearing. 

The rotor is driven at uniform speed by a 
one-half horse power induction motor. The speed 
may be varied continuously from 360 r.p.m. to 
3600 r.p.m. by a DT-} Transitorq unit which is 
attached directly to the drive shaft of the motor. 
Lower speeds are made available by a 3-1 
reduction pulley. The power is transmitted from 
the Transitorq unit to the rotor by a flexible 
shaft whose length is 6’ and core diameter is 2”. 
A stroboscope is mounted at the Transitorq end 
of the flexible shaft. Its accuracy in determining 
speeds of rotation, as checked against a tacho- 
meter, is 0.5 percent. 

The optical system consists of light source, 
filter, slit, collimating lens, polarizing Nicol, 
Kundt cell, mica quarter-wave plate compen- 
sator, mica half-wave plate Chauvin half-shadow, 
analyzing Nicol, and telescope. 

The light source is Type H-3 85-watt General 
Electric Company mercury vapor lamp. The 
mercury green line (5461A) is isolated by Corning 
glass filters Nos. 351, 430, and 512. The filters 
are mounted on the stand supporting the slit 
which is usually set at 0.04-0.06 cm. The lens 
is placed at its focal length, 15 cm, from the 
slit to obtain a parallel beam. 

The polarizing and analyzing Nicols of the 
polarimeter may be locked together in the 
crossed position and rotated as a unit for the 
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purpose of determining the vortical angle. The 
coupling is made by two identical brass drums 
which operate on a drill rod steel shaft. These 
drums are connected by piano wire to two larger 
identical drums which rotate with the Nicol 
prisms. Turning the steel shaft causes the two 
Nicols to rotate without changing appreciably 
their orientation over a rotation of 90°. 

Since the vortical angle is referred to the 
direction of flow, and the latter is horizontal 
when the Kundt cell is operated with its windows 
at the top, the horizontal position of either 
polarizer or analyzer must be found. This 
calibration was made by suspending a glass plate 
vertically in front of the analyzer. Then the 
planes of the glass and light source were adjusted 
to reflect vertically plane polarized light into 
the analyzer, which was rotated to extinction. 
This experiment determined the reading on the 
polarimeter scale at which the principal plane of 
the analyzing Nicol is horizontal. 

The assembled cell is filled with solution 
through the filling tube, by a capillary attached 
to a 40-cc separatory funnel which serves as 
reservoir for the solution during an experiment. 
Bubbles are flushed out through two small 
holes drilled through the window holders. These 
holes are closed with wooden plugs when the 
cell is in operation. When filled, the cell is 
clamped horizontally on the optical bench and 
connections are made to the water lines and 
flexible shaft. The Transitorq is set at any 
desired speed. The vortical angle y is determined 
by removing the compensator and rotating the 
locked crossed Nicols from the reference hori- 
zontal position to the nearest minimum. To 
measure the amount of birefringence, the rotor 
is stopped, and the crossed Nicols are rotated 45° 
from the minimum just found. The compensator 
is inserted and rotated to extinction. The polar- 
izer is locked in position and the analyzer un- 
coupled. The rotor is again started and the 
analyzer alone is rotated to extinction. Thus, 
rotation of the plane of the polarized light in 
degrees is measured on the polarimeter scale. 
It is the angle Z of Eq. (1). 

The accuracy of the determination of the 
vortical angle depends upon the amount of 
birefringence involved, siice this determines the 
contrast between maxima and minima observed 
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as the crossed Nicols are rotated. When the 
angle Z is less than 0.5°, the maxima and minima 
are too broad to distinguish. For values of Z in 
the vicinity of 1.0°, y may be determined with a 
probable error of 5° or less. The angle Z itself 
may be determined with an accuracy of +0.05° 
under the best conditions. Eye fatigue, or any 
factors which cut down the light intensity, will 
reduce this accuracy. 


4. PREPARATION OF MATERIALS 


A series of polymers of a-methyl methacrylate 
was prepared by heating the monomer in solution 
according to a method suggested by Strain.!° 
The reaction was carried out under a number of 
different conditions to obtain products of various 
degrees of polymerization. The relative degree of 
polymerization was roughly estimated by meas- 
uring the viscosity of 0.1 percent solutions of the 
products in dried dioxane. Of nine polymers pre- 
pared, five were used for the double refraction ex- 
periments. These will be designated in order of 
increasing molecular weight as A, B, C, D, and E. 

Since it is essential that the samples be as 
homogeneous as possible, various techniques 
were employed to reduce the polydispersity. 
The degree of polymerization was found to 
increase during the course of a reaction. Hence, 
one useful technique involved the removal of 
the product at various stages of the reaction. 
This was possible in the case of the last four 
polymers because the solvent used, aqueous 
methanol, precipitates the polymer as it is 
formed. 


Polymer A 


Ten cc of monomer and 40 cc of commercial 
dioxane were heated at 65°+1°C under reflux. 
Ten cc of monomer was added at a rate of 1 cc 
per five hours. This rate corresponded to the 
initial rate of the reaction as determined by 
preliminary experiments. By keeping the mono- 
mer concentration constant in this way it was 
hoped to decrease the polydispersity of the 
product. When all the monomer had been added, 
the solution was diluted with a large excess of 
acetone, and the polymer was precipitated with 
water. It was then washed with water and dried 
at 85°C. 


10 Strain, Ind. Eng. Chem. 30, 345 (1938). 
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Polymers B and C 


Fifty cc of monomer and 200 cc of 50 percent 
by volume aqueous methanol, together with 1 
percent of benzoyl peroxide figured on the basis 
of the monomer were heated under reflux at 
65°C. After an induction period of about two 
hours, precipitation commenced. At the end of 
three hours the precipitate was removed, washed 
with methanol and water, and dried. This frac- 
tion represented a portion of the total yield 
included between 0 percent and 49 percent. It is 
designated as polymer B. 

The reaction was repeated, the cloudy super- 
natant liquid being decanted into a clean re- 
action flask at the end of 2} hours. At the end 
of an additional 14-hour period, more precipitate 
had formed in the fresh flask. This material, 
designated as polymer C, was washed with 
methanol, water, and dried. It represented a 
fraction of the total yield included between 50 
percent and 80 percent. 


Polymers D and E 


The conditions for the preparation of these 
polymers were the same as for polymers B and C, 
except that the catalyst, benzoyl peroxide, was 
omitted. Precipitation began at the end of 50 
hours, and by 64 hours the reaction had gone 
to completion. The total product was separated 
into two parts, a hard central lump, and a soft 
material which had accumulated on the walls of 
the flask. The latter fraction, designated as 
polymer D, was washed with water and dried. 
The aggregated material was further fractionated 
by dissolving it in acetone and diluting the 
solution with two volumes of hot ethanol. As 
the solution cooled to room temperature, the 
polymer precipitated. After four hours the pre- 
cipitate which had settled was removed, washed, 
and dried. This fraction is designated as poly- 
mer E. 

The dioxane which was used as solvent in the 
double refraction measurements was carefully 
purified according to a procedure suggested by 
Weissberger and Proskauer."' The material gen- 
erally showed upon distillation a boiling range 
of less than 0.3°. Ten samples gave d,25= 1.0273, 
and 725=0.01186 as average values. 


1! Weissberger and Proskauer, Organic Solvents (Claren- 
don Press, Oxford, 1935). 
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DOUBLE REFRACTION 


The viscosities of 0.1 percent solutions of 
polymers A, B, C, D and E in dioxane at 25°C 
were 0.0125, 0.0140, 0.0147, 0.0162 and 0.0191, in 
order. They are here expressed in absolute units. 


5. EXPERIMENTAL RESULTS 


Observations of vortical angle and birefringence 
were made at 25° for dioxane solutions of the 
five methyl methacrylate polymers over a range 
of concentration from 1 gram per liter to 50 
grams per liter, and over a range of velocity 
gradient from 4000 sec.— to 15,000 sec.—. 

In all cases the magnitude of the birefringence 
was found to be proportional to the velocity 
gradient. This law has been shown to hold for 
a large number of systems (polystyrene, nitro- 
cellulose,!? w-hydroxy decanoic acid polyesters," 
tobacco mosaic virus protein,’ and a large 
number of pure liquids.'® Such behavior is pre- 
dicted by Boeder in his solution of Eq. (3), 
provided the ratio G/@ is small. 

From this fact it follows that the birefringence 
S previously defined is a constant for any given 
solution. In Fig. 1, S is plotted for each of the 
five polymers as a function of concentration. 
The graph shows that in dilute solutions S$ 
varies linearly with concentration, and that all 
the lines extrapolate to the point marked Sp at 
infinite dilution. This point should represent the 
magnitude of the solvent birefringence. Actually, 
So for dioxane was found to be 0.1310- by 
averaging the results of 57 measurements taken 
at a velocity gradient of 6200 sec.—'. This value 
agrees within experimental error with that ob- 
tained by extrapolating S to zero concentration. 
It is concluded, therefore, that the specific 
birefringence D is a constant characteristic of 
each polymer in a given solvent. 

Largely as a matter of convenience in the 
evaluation of the data, the plots of Fig. 1 have 
been drawn as broken lines. The slopes of the 
S us. ¢ curves are constant in regions of low con- 
centration and give the specific birefringence D 
for the several polymers, after correction has 
been made for the solvent birefringence (Eq. (2)). 


aad and Gross, Zeits. f. physik. Chemie A165, 161 





(1 
8 Quimby, dissertation, University of Wisconsin (1937). 
Mehl, Cold Spring Harbor Symposium on Quantitative 

Biology, Vol. 6 (1938), p. 218. 

** Sadron, J. de phys. et rad. 7, 263 (1936). 
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OF FLOW 771 
The specific birefringence ceases to be constant 
in the more concentrated solutions and it in- 
creases with concentration. (Additional data for 
polymers A and E are not shown in the figure 
for lack of space, but are indicated by the direc- 
tion of the arrows.) 

Preliminary to any measurement of B, the 
vortical angle y must be determined. Table II 
gives the average values of the vortical angles for 
solutions of the five polymers, measured over a 
range of velocity gradients. The data in italics 
are the G/@ ratios corresponding to each vortical 
angle, as interpolated from Table I. 

One would expect, from Boeder’s theory, that 
the vortical angle should be independent of 
concentration. Actually it shows a marked con- 
centration dependence, a phenomenon which was 
also observed by von Muralt and Edsall'® in 
working with muscle globulin. However, it has 
just been demonstrated that the specific bire- 
fringence varies in solutions of appreciable con- 
centration. Therefore, by Eqs. (5) or (7), @ should 
vary and hence also y. With this consideration in 
mind we have selected for the average only 
those values of Y which have been measured in 
the concentration region where D is constant. 
Unfortunately, in this concentration region where 
one would expect y to remain constant, experi- 
mental observations are less accurate, and its 
exact behavior cannot be determined with the 
present apparatus. One factor which might well 
cause variation in the vortical angle even below 
the critical concentration is the additive effect 
of the vortical angle of a slightly birefringent 
solvent, which is always 45°. Sadron'? has shown 
that for such a two-component system 


tan 2Wons — tan 2y1= (Bo/B,) sec 21 (8) 


where the subscript 1 refers to the solute and 
subscript 0 to the solvent. If Bo is less than 75 
the observed birefringence, then the effect of the 
solvent on the vortical angle is within the 
experimental error. 

The rotary diffusion constants of the five 
polymers in radians per second were calculated 
by dividing the G/@ ratios (Table II) into the 
corresponding values of the velocity gradient 


16 yon Muralt and Edsall, J. Biol. Chem. 89, 315, 351 
(1930). 
17 Sadron, J. de phys. et rad. 9, 381, 384 (1938). 
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Fic. 1. Birefringence of methyl methacrylate polymers in 
dioxane solution, per unit velocity gradient, per unit 
viscosity of solvent as a function of concentration. 


and averaging the results. In Table III they have 
been reduced to unit viscosity by multiplying 
by 0.0119, the absolute viscosity of the solvent. 

In the course of assembling the birefringence 
data, the viscosity of each solution was measured 
at 25°C. The specific viscosity increment at 
infinite dilution was calculated for each polymer 
by extrapolation. In Table III, the symbol C is 
used to indicate that the concentration is ex- 
pressed in volume fraction rather than in grams 
per liter. The density of the polymers was taken 
as 1.18 in making the conversion. 


5. Discussion 


The viscosity correction 


The specific birefringence D has been defined 
by Eq. (2) with reference to the viscosity of 
the solvent. Signer and Gross!* have reported 
that, for low molecular weight polystyrenes, the 
specific birefringence does not remain constant 
even in the most dilute solutions studied, unless 
it is referred to the viscosity of the solution. 
Mehl" discusses this point in connection with 
the double refraction of flow of tobacco mosaic 
virus protein. His own data are not extensive 
enough to permit a choice between the viscosity 
of solvent and solution in calculating D. In 
reviewing the experiments of von Muralt and 
Edsall on muscle globulin, however, he finds 
that the viscosity of the solution gives the more 
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constant values for D. Since the equations in- 
volved in the theory of double refraction of flow 
and of the rotary diffusion constant refer to the 
viscosity of the solvent, empirical relationships 
involving the viscosity of the solution are diffi- 
cult to interpret. It is probable that the correc- 
tion can be avoided by working at lower con- 
centrations. 

Figure 1 shows that for dioxane solutions of 
methyl methacrylate polymers, D is constant and 
independent of the viscosity of the solution over 
a considerable range of molecular sizes. The 
failure of Signer and Gross to obtain a similar re- 
sult may be attributed largely to their neglect of 
the birefringence of the solvent, cyclohexanone. 
Recalculation of their data for two polystyrenes 
of molecular weight 30,000 and 80,000 indicates 
that curves of the type shown in Fig. 1 may be 
obtained for these macromolecules. In the case 
of the muscle globulin experiments, it is quite 
possible that the concentrations used, 3 grams 
per liter to 15 grams per liter, belong in the region 
where D varies with concentration. The variation 
of y with concentration found in these experi- 
ments indicates that such is actually the case. 


The relationship between D and 4 


In Fig. 2, the reciprocal of the rotary diffusion 
constant, reduced to unit viscosity, is plotted as 
a function of the specific birefringence. The 
points, except for that given by polymer C fall 
fairly well on a straight line. The data can be 
compared with the requirements of the theo- 
retical equations (5) and (7). Equation (5) de- 
TABLE II. Vortical angles of dioxane solutions of methyl 


methacrylate polymers. (Corresponding values of 
G/@ are given in italics.) 








GRADIENT (SEC.~) 





PoLYMER 4150 6200 8300 10,300 12,500 15,500 
A 37.0 35.0 
— — —_ _ 1.15 1.50 
B 31.5 28.0 26.5 
— — — 2.15 2.90 3.30 
Cc 34.0 30.5 
— — — ~_ 1.70 = 2.35 
D 29.5 23.0 25.0 22.5 17.0 
_— 2.55 4.50 380 4.75 7.95 
E 254 206 22.5 20.7 17.0 
3.70 5.75 4.75 5.75 7.95 = 
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mands that the data fall on a straight line which 
passes through the origin. The relationship of 
Eq. (7) is not so obvious. Sadron has shown that 
for the polystyrenes, A’. increases with molecular 
weight. Assuming this to hold true for methyl 
methacrylate polymers, and remembering that 
6) decreases with increasing molecular weight, 
one sees that the variation in D due to the term 
in 1/6) is of higher order than the variation due 
to A’/@.A’.. Hence, Eq. (7) may be written in 
the form, 

D=k,[1/00—ke], (9) 


where k; is a positive empirical constant, and 
where the positive quantity, k2=A’/@A’., is 
assumed as a first approximation to be constant. 
The curve of Fig. 2 is consistent with the require- 
ments of Eq. (9). 

While the data for the vortical angle are 
probably not accurate enough to distinguish 
between Eqs. (5) and (7), the latter has a more 
sound theoretical basis, and the rotary diffusion 
constant will be considered to obey the approxi- 
mate law of Eq. (9). On this assumption one 
may refer the rotary diffusion constants of the 
five polymers to their specific birefringences, 
rather than to the less accurately measurable 
vortical angles. Values of 4, corrected on this 
basis, are given in Table III. 


Mass, shape, and structure of the methyl 
methacrylate polymers 
From a priori considerations of the structure 
of methyl methacrylate polymers, one might 
predict that the macromolecule will behave as a 
rigid rod in solution. If one assumes the structure 
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TaBLe III. Birefringence and molecular shape data for 
methyl methacrylate polymers (25°). 








d 
Po.y- oie nsp/C 





60 60* Ss 
MER D X10 Oss. Corr. (A) M Ratio) (C0) 
A 0.09 127 149 890 34,000 89 75 
B- 0.29 55 64 1,170 45,000 117 210 
C 0.42 84 48 1,300 50,000 130 270 
D_ 0.69 27.5 32 1,500 58,000 150 410 
E 1.32 17.5 18 1,880 72,000 188 640 








* Rotary diffusion constants are measured in radians per second at 
unit viscosity. 
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Fic. 2. Reciprocal of rotary diffusion constant as a 
function of specific birefringence for polymers A, B, C, D 
and E (cf. Table III). 


then the strongly polar ester groups occur on 
alternate carbon atoms. At such a small distance, 
the repulsion between ester groups will be great, 
and will minimize any tendency towards “‘kink- 
ing,’ or rotation about the carbon-carbon bonds 
in the chain. This structure may be contrasted 
with that of w-hydroxydecanoic acid polyester, 
in which a chain of nine single bond carbon to 
carbon linkages separates each polar group. 

A Fischer-Hirschfelder molecular model of 
polymethyl methacrylate was constructed to 
correspond with the above structure. From an 
inspection of this model it was concluded that 
the macromolecule acts as a rigid rod, since there 
is little opportunity for kinking in the model. 

Since Eq. (4) applies specifically to a long, 
rigid macromolecule, it can be used to calculate 
the lengths of the five polymers from their 
rotary diffusion constants. The corrected values 
of the latter were used in the calculation. The 
results are given in Table III. 

The molecular model indicates that the poly- 
methyl methacrylate macromolecule has a mean 
diameter of 10A. This will depend upon the 
particular structure assigned to the “‘backbone”’ 
carbon chain of the polymer. In our model, this 
chain consisted of a single row of carbon atoms. 
A chain consisting of a double row of carbon 
atoms is also sterically permissible, and this 
model would indicate a larger diameter. The 
values of the ratio of length to width s/d are 
listed in Table III for d=10A. 
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One may calculate the molecular weights of the 
polymers on the basis of a single row ‘“‘backbone”’ 
chain by multiplying the values of s from Table 
III by 100/2.6, where 100 is the molecular 
weight of the monomeric unit and 2.6A is the 
length which the unit occupies in the chain. 
The molecular weights so calculated, also given 
in Table III, are minimum values, and again 
they depend upon the structure assumed for the 
carbon chain. The relative values of the mo- 
lecular weight should be, of course, independent 
of this assumption. 


Correlation with the specific viscosity increment 


In Fig. 3 are plotted the specific viscosity 
increments, 7s»/C (extrapolated to infinite dilu- 
tion) as a function of the square of the length to 
width ratio of the polymers. The dotted lines 
are obtained from theoretical equations (Kuhn, 
Guth-Eisenschitz). The dependence of the specific 
viscosity increment on the square of the length 
to width ratio is well substantiated. The slope 
of the experimental line is sensitive to the value 
assumed for the diameter of the macromolecule. 
If d2=10A as assumed, the slope does not differ 
greatly from that demanded by the Guth-Eisen- 
schitz equation. If d is assumed to be 16A, the 
slope corresponds to that given by the Kuhn 
equation. The most important discrepancy be- 
tween theory and experiment is that the experi- 
mental line does not pass through the origin. 
The reason for this disagreement is not clear, 
although it is possible the failure to extrapolate 
specific viscosities to zero shear is involved. 
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During the course of the experimental work it 
was assumed the flow is Newtonian because of 
the low concentrations of the solutions and the 
low gradients used. 


Spontaneous orientation 


It was pointed out in the previous section 
that the specific birefringence of the methyl 
methacrylate polymer solutions increases after 
a region of higher concentration is reached. On 
the basis of the orientation theory, this increase 
may be interpreted as being due to the action of 
orienting forces within the solution which are 
negligible at low concentrations, but which be- 
come operative when the solute concentration is 
increased. 

The exact nature of these forces which produce 
spontaneous orientation is not definitely known. 
They may consist partly of long range van der 
Waals forces of the type described by London.'* 
The macromolecules may tend to orient simply 
because this configuration corresponds to a 
minimum volume, and hence to a minimum 
energy for the system. The latter effect should 
disappear at low concentrations where no close 
packing is involved. 

The data obtained for methyl methacrylate 
polymers may give positive correlation between 
a concentration at which spontaneous orientation 
sets in and the size of the molecule. If it is not 
unreasonable to interpret the data as has been 
done in Fig. 1, where all the curves turn upwards 
at the same value of S, the critical concentration 
is approximately proportional to the rotary 
diffusion constant, or to the inverse cube of the 
molecular length. 
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18 London, Trans. Faraday Soc. 33, 8 (1937). 
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The Exchange Reaction between Gaseous and 
Combined Nitrogen 
GeorGEs G. Joris 


Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 
September 9, 1941 


HE exchange reaction between gaseous and combined 
nitrogen recently reported in a letter to this journal 
by Yoshio Nishina, Takeo Iimori, Hideo Kubo and Hiromi 
Nakayama! has been reinvestigated using stable heavy 
nitrogen isotopes N*, and N?%. The composition 
of the initial gaseous mixture was N?*%.: N#%, : N88, 
1000 : 204 : 100. Aqueous solutions of NaNO, NaNO; and 
NH-OH-HC1 at a concentration of 1/20 mole/l were used 
since they were reported to undergo the most rapid ex- 
change between gaseous nitrogen and combined nitrogen. 
No exchange was found as shown by the ratios N?%2/N#°, 
with the shaking time indicated in Table I. 











TABLE I. 

SUBSTANCE TESTED 0 HOUR 6 HOURS 30 HOURS 360 HOURS 
NaNOs 2.04 2.05 2.05 2.04 
NaNOz 2.04 2.05 2.04 2.05 
NH:OH -HCl 2.04 2.03 2.05 2.03 








The results recorded with radioactive nitrogen must 
therefore be ascribed to difficulties arising in the technique 
of operation with the radioactive isotope. The present 
results are more reasonable, judged from the chemical 
standpoint and the stability of the ions containing nitrogen. 


1 Y. Nishina, T. Iimori, H. Kubo and H. Nakayama, J. Chem. Phys. 
9, 571 (1941). 





The Infra-Red Absorption and Reflection 
Spectra of KHF, and of KDF; 


J. A. A. KETELAAR 


Laboratory for Inorganic and Physical Chemistry, 
University, Leyden, The Netherlands 


September 5, 1941 


HESE compounds were investigated with a prism 
spectrometer in the quartz and rocksalt regions from 
1-164. When the investigation of ordinary potassium bi- 
fluoride was about finished a note came to my knowledge 
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from Messrs. Buswell, Maycock and Rodebush.' These 
authors only investigated the quartz region and here our 
results are incomplete accordance. However my much more 
extensive investigations have led to conclusions essentially 
different from theirs. 

The maxima at 1222 cm™ and 1450 cm™ (in reflection 
1238, 1490 cm~) are both fundamental bands as indicated 
by the extreme strong absorption (> 90 percent for a layer 
of 84) and by the high reflection percentages (55 percent 
and 60 percent, respectively) for both bands, thus excluding 
the explanation as a combination frequency for the latter 
band. The three other pairs of bands of decreasing intensity 
below 3500 cm~ are due to combination frequencies of 
these fundamentals with 1, 2, and 3 times a frequency of 
560-600 cm~!. Of these pairs also the higher frequency 
component is the strongest band. The bands above this 
limit are due to the second harmonics? uf the fundamentals 
at 3730 cm™ (double) and at 5099 cm™, respectively, and 























Fic. 1. Transmission and reflection curves for KHF2 
in the rocksalt region. 


to their combinations with the same frequency of 560 cm™. 
The doubling of the band at 3730 cm™, also observed by 
the other authors, is due to a minor complication, the 
splitting of an higher harmonic due to higher terms in the 
potential function. The spectrum of KDF: is explained in 
just the same way with a shift of the fundamental bands 
by about v2 to 891 cm™ and 1046 cm™ (in reflection). 
The inactive frequency at 560 cm™ is found to be inde- 
pendent of the mass of the hydrogen atom. Thus a com- 
plete analysis of all bands in both spectra is possible with a 
single exception of minor importance. 

The spectrum observed is due to the isolated HF» group 
as indicated by the identity between the spectra of KHF»2 
and RbHFs», save for a shift of a few cm™ to lower values 
for the latter compound. Moreover specimens of KDF:2 con- 
taining 10 percent H give a mere superposition without any 
shift of the spectra of the pure compounds, thus proving 
the absence of any appreciable coupling between the 
groups. The fundamental double band is to be associated 
with the asymmetric valency oscillation of a linear three 
atomic model (‘‘hydrogen frequency”), whereas the fre- 
quency at 560 cm is the symmetric one (‘‘fluorine fre- 
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quency”). The doubling of the active frequency can be 
shown to be only due to the existence of a double minimum 
in the potential carve for the hydrogen atom. A potential 
curve constructed according to the method of Dennison 
and Uhlenbeck® from the experimental data gives 2570 
cm for the height of the potential hump with a separation 
of both levels of 25 cm and 185 cm™ for the ground and 
for the first excited level, respectively, the latter being 
already very near the top of the hill. The distance between 
both minima is found as only 0.70 angstrom units. 

These results prove for the first time the existence of 
protonic resonance in the case of a hydrogen bond. On the 
other hand the small doubling of the ground level proves 
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definitely that even here this resonance does not contribute 
to any extent to the energy of the bond formation (hy 
=0.04 kcal./mole) which consequently must be considered 
to be entirely of electrostatic origin, due to the partial 
ionic character of the hydrogen atom. The structure of the 
HF: group is thus symmetric, resonating between F HF 
and FH F. 

Experimental details are given in a forthcoming publica- 
tion in the Rec. trav. chim. 

1A. M. Buswell, R. L. Maycock and W. H. Rodebush, J. Chem. 
Phys. 8, 362 (1940). 


2 The first harmonics are absent due to a selection rule valid for the 


crystal. 
3D. M. Dennison and G. E. Uhlenbeck, Phys. Rev. 41, 313 (1932). 








